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Abstract 

The  diffraction  of  a  wave  by  an  aperture  of  any  shape  in  a  thin  screen 
is  treated  by  a  new  method.  This  method  is  called  "  the  geometrical  theory  of 
diffraction",  because  it  is  an  extension  of  geometrical  optics  which  accounts 
for  diffraction.  In  this  method  new  rays  -  called  diffracted  rays  -  are  intro- 
duced. They  are  produced  when  an  incident  ray  hits  the  edge  of  the  aperture, 
and  they  satisfy  the  "  law  of  diffraction".   A  field  is  associated  with  each 
ray  in  a  quantitative  way,  by  means  of  the  optical  principles  of  phase  variation 
and  energy  conservation.  In  addition  "  diffraction  coefficients  n  are  introduced 
to  relate  the  field  on  a  diffracted  ray  to  that  on  the  corresponding  incident  ray. 

Ey  this  method  a  simple  formula  is  obtained  for  the  field  diffracted 
by  any  aperture.  By  means  of  this  formula  the  field  in  the  aperture  or  the  far 
field  diffraction  pattern  can  be  found.  In  addition  the  transmission  cross 
section  of  the  aperture  can  be  determined.  Explicit  formulas  and  numerical  re- 
sults are  given  for  slits  and  circular  apertures.  The  accuracy  of  the  results 
increases  as  the  wavelength  decreases,  but  they  are  even  useful  for  wavelengths 
as  large  as  the  aperture  dimensions. 
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1.  Introduction 

The  present  treatment  of  diffraction  of  a  wave  by  an  aperture  in  a  thin 
screen  is  based  upon  a  new  theory  of  wave  propagation.  This  theory  is  called 
•the  geometrical  tneory  of  diffraction' ^  J,  because  it  is  an  extension  of 
geometrical  optics  which  describes  diffraction  phenomena.  It  is  like  geometri- 
cal optics  inasmuch  as  it  postulates  that  fields  propagate  along  certain  lines 
or  curves  called  rays.  However,  it  goes  beyond  geometrical  optics,  by  intro- 
ducing various  new  kinds  of  rays,  called  'diffracted  rays',  in  addition  to  the 
usual  optical  ones.  It  also  assigns  a  field  value  to  each  point  on  a  ray  in  a 
quantitative  way  and  postulates  that  the  total  field  at  a  point  is  the  sum  of  the 
fields  on  all  rays  through  that  point.  Thus  the  new  theory  not  only  provides  a 
qualitative  description  of  diffraction  in  terms  of  the  diffracted  rays,  but  it 
permits  a  quantitative  determination  of  the  field  as  well. 

By  using  this  theory  we  have  obtained  a  relatively  simple  formula  for 
the  field  diffracted  by  an  aperture  of  any  shape.  This  formula  takes  into  account 
the  nature  of  the  field  under  consideration,  i.e.  whether  it  is  electromagnetic, 
acoustic  or  of  some  other  type,  and  the  nature  of  the  screen.  It  enables  us  to 
compute  the  field  in  the  aperture  as  well  as  in  front  of  and  behind  the  screen. 
In  particular  it  yields  the  diffraction  pattern  and  the  transmission  cross  section 
of  the  aperture.  The  formula  also  applies  to  diffraction  by  a  flat  thin  obstacle 
of  any  shape,  and  thus  permits  a  verification  of  Babinet's  principle.  It  also 
accounts  for  the  bright  lines  and  bright  spots  which  have  been  observed  in  the 
shadows  of  such  obstacles. 

Because  our  theory  is  related  to  geometrical  optics,  our  results  should 
be  valid  when  the  wavelength  X  is  small  compared  to  the  aperture  dimensions.  We 
have  verified  this  expectation  by  comparing  our  results  with  certain  previously 
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known  exact  results.  In  doing  so  we  have  also  found  that  our  results  are 
accurate  enough  to  be  of  use  even  for  wavelengths  as  large  as  the  aperture 
dimensions.  Of  course  accuracy  increases  as  the  wavelength  decreases. 

Previously  we  have  used  the  'geometrical  theory  of  diffraction'  to 

[2l 
treat  diffraction  by  a  cylinder  of  any  cross  section1-  J ,  diffraction  by  a  smooth 

three-dimensional  object  of  any  shape'-  J,  and  other  diffraction  problems'-  K     In 
these  cases  we  also  found  excellent  agreement  between  our  results  and  previously 
known  exact  results.  In  fact,  whenever  comparison  was  possible,  our  results 
agreed  perfectly  with  the  leading  terms  in  the  asymptotic  expansions  of  the 
exact  solutions  of  diffraction  problems,  the  expansions  being  with  respect  to 
k  =  2ti A>  as  k  ->  co .  This  agreement  has  confirmed  our  belief,  which  is  based 
upon  various  mathematical  considerations,  that  our  theory  always  yields  the  lead- 
ing terms  in  the  asymptotic  expansion  of  the  exact  solution.  Since  such  expansions 
are  often  accurate  even  for  small  values  of  k  (large  values  of  X),  the  accuracy  of 
our  results  for  large  wavelengths  is  not  surprising. 

Certain  aspects  of  our  theory  are,  like  geometrical  optics,  the  same 
for  many  kinds  of  fields.  This  is  true  of  the  rules  governing  the  rays,  both 
ordinary  and  diffracted,  and  the  principles  determining  the  variation  of  phase 
and  amplitude  along  a  ray.  Therefore  these  parts  of  the  theory  can  be  presented 
in  a  general  manner,  without  reference  to  the  nature  of  any  particular  field. 
Gn  the  other  hand  there  are  other  aspects  of  the  theory  which  depend  essentially 
upon  the  nature  of  the  field.  This  is  the  case  with  coefficients  of  reflection 
and  transmission,  and  with  the  diffraction  coefficients  which  we  will  introduce. 
Consequently  these  coefficients  must  be  considered  separately  for  different  kinds 
of  fields.  Decause  of  this  we  will  first  present  the  theory  for  a  scalar  field 
(e.g.  acoustic  pressure)  and  later  show  what  differences  occur  in  the  treatment 
of  electromagnetic  fields. 
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We  will  formulate  our  problem  as  follows:  A  thin  screen  containing 
an  aperture  lies  in  the  plane  x  =  0 .  The  rest  of  space  is  f illed  with  an 
homogeneous  medium  with  propagation  speed  c.  A  given  wave,  with  wavelength  \ 
and  propagation  constant  k  «  •=—  >   is  incident  upon  the  screen  from  the  left.  Vfe 
wish  to  find  the  field  which  results  from  the  interaction  between  the  incident 
wave  and  the  screen. 

The  incident  wave  can  be  characterized  in  either  of  two  ways.  The  first 
characterization  consists  in  giving  the  set  of  incident  rays  and  the  field  value 
(i.e.  phase  and  amplitude)  on  each  ray.  The  second  alternative  consists  in 
giving  the  value  of  the  incident  field  u.   at  each  point  in  the  form 

-Lric 


(1)         Hf„>.y.»)   c  A(x,y,«)  e1"^'^ 


inc 


From  the  second  form  (l)  the  rays  can  be  determined  as  the  orthogonal  trajectories 
of  the  wavefronts  Y(x,y,z)  ■  constant.  Conversely,  when  the  rays  and  the  phase 
on  each  ray  are  known,  Y(x,y,z)  can  be  found.  Thus  the  two  characterizations  are 
equivalent.  Therefore  we  will  assume  that  u.   (x,y,z)  is  given  in  the  form  (l) 
and  that  the  incident  rays  are  also  known. 

In  Section  2  we  will  determine  the  rays — especially  the  diffracted  ones-- 
which  are  produced  when  the  incident  rays  hit  the  screen.  In  Section  3  we  will 
determine  how  the  field  varies  along  a  ray.  In  S3Ction  h  we  will  determine  the 
initial  value  of  the  field  on  each  ray.  Then  in  Section  5  we  will  obtain  the 
formula  for  the  field  at  any  point  by  adding  together  the  fields  on  all  rays 
through  the  point.  In  subsequent  sections  the  results  will  be  applied  to 
diffraction  by  a  half-plane,  a  slit,  a  circular  aperture,  a  convex  aperture, 
and  a  grating  of  slits. 

Most  of  the  previous  work  on  diffraction  by  an  aperture  has  been  carefully 
discussed  by  C.J.  Bouwkamp'-  ■* , 
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2.  Ray  tracing 

According  to  geometrical  optics,  the  only  rays  which  occur  in  the 
present  problem  are  the  incident  rays  and  the  reflected  rays  to  which  they  give  rise 
when  they  hit  the  screen.  All  of  them  are  straight  lines,  and  the  direction  of 
each  reflected  ray  is  determined  by  the  law  of  reflection. 

V/e  assume,  however,  that  there  are  additional  rays  which  are  produced  when 
incident  rays  hit  the  edge  of  the  screen.  We  call  them  'diffracted  rays'  and  assume 
that  they  are  straight  lines  and  that  they  satisfy  the  'law  of  edge  diffraction. ' 
This  law  states  that  a  diffracted  ray  and  the  corresponding  incident  ray  make 
equal  angles  with  the  edge  and  lie  on  opposite  sides  of  the  plane  normal  to  the 
edge  at  the  point  of  diffraction.  From  this  law  it  follows  that  there  are 
infinitely  many  diffracted  rays  produced  by  one  incident  ray.  These  diffracted 
rays  lie  on  the  surface  of  a  cone  having  as  axis  the  tangent  to  the  edge  at  the 
point  of  diffraction  (see  Figure  l).  The  half-angle  of  the  cone  is  just  the 
angle  between  the  incident  ray  and  the  edge.  When  the  incident  ray  is  perpendicular 
to  the  edge,  the  half-angle  is  n/Z   and  the  cone  of  diffracted  rays  becomes  a  plane 

perpendicular  to  the  edge. 

The  'law  of  edge  diffraction'  and  the  straightness  of  the  diffracted  rays 
are  consequences  of  an  extended  form  of  Fermat's  principle W.  The  special  form 
of  this  principle  which  suffices  in  the  present  problem  is  the  following:  A 
singly  diffracted  ray  connecting  two  points  P  and  Q  is  a  curve  whose  length  is 
stationary  among  all  curves  connecting  these  two  points  and  having  one  point  on 
the  edge  of  the  screen.  From  this  principle  it  follows  that  a  diffracted  ray 
consists  of  two  straight  line  segments,  one  from  P  and  one  from  0  to  some  point 
on  the  edge.  Furthermore  these  lines  must  make  equal  angles  with  the  edge  and 
lie  on  opposite  sides  of  the  plane  normal  to  the  edge  at  the  point  of  diffraction. 
This  result  is  the  law  of  edge  diffraction  stated  above. 
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Figure  la. 

The  cone  of  diffracted  rays  produced 
by  an  incident  ray  which  hits  the  edge 
of  a  thin  screen. 


Figure  lb. 

The  plane  of  diffracted  rays 
produced  by  a  ray  normally 
incident  on  the  edge  of  a  thin 
screen. 
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We  also  assume  that  a  diffracted  ray  which  hits  the  edge  can  again 
produce  diffracted  rays  in  accordance  with  the'law  of  edge  diffraction1.  Thus 
doubly,  triply  and  multiply  diffracted  rays  can  occur.  All  of  these  rays  also 
follow  from  the  extended  Fermat's  principle  in  which  curves  with  more  than  one 
point  on  the  edge  are  considered. 

Aa  an  example  of  the  preceding  considerations  Figure  2  shows  the  rays 
produced  by  a  plane  wave  normally  incident  upon  a  thin  screen  containing  a  slit. 
Since  each  incident  ray  which  hits  an  edge  of  the  screen  is  perpendicular  to  this 
edge,  the  resulting  diffracted  rays  lie  in  a  plane  normal  to  the  edge.  Through 
each  point  behind  the  screen  there  pass  two  singly  diffracted  rays  (one  from  each 
edge),  two  doubly  diffracted  rays,  etc.  In  addition  an  incident  ray  passes 
through  each  point  not  in  the  shadow. 

So  far  we  have  not  considered  any  ray  which  hits  a  corner  or  vertex  of 
the  edge.  We  assume  that  such  rays  also  give  rise  to  diffracted  rays  which  are 
straight  lines  satisfying  the  'law  of  vertex  diffraction. '  This  law  states  that 
a  ray  diffracted  from  a  vertex  may  have  any  direction,  without  any  relation  to  the 
direction  of  the  incident  ray.  As  a  consequence  of  this  law  it  follows  that  a 
single  ray  which  hits  a  vertex  produces  diffracted  rays  which  leave  the  vertex  in 
all  directions. 

The  straightness  of  the  vertex  diffracted  rays  and  the  law  of  vertex 
diffraction  also  follow  from  the  extended  form  of  Fermat's  principle"-  -I .  The 
special  form  of  this  principle  which  characterizes  the  vertex  diffracted  rays  is 
this:  A  ray  connecting  two  points  P  and  C  and  singly  diffracted  from  the  vertex  V 
is  a  curve  whose  length  is  stationary  among  all  curves  connecting  these  two  points 
and  passing  through  V.  This  principle  implies  that  such  a  ray  consists  of  two 
straight  line  segments,  one  from  P  to  V  and  one  from  V  to  C.  The  directions  of  these 
two  lines  are  unrelated. 

A  ray  diffracted  from  a  vertex  may  hit  another  vertex  or  an  ede;e  and 
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Figure  2 

The  diffracted  rays  produced  by  a  plane  wave  normally 
incident  upon  a  slit  in  a  thin  screen.  The  two  incident 
rays  which  hit  the  slit  edges  are  shown  along  with  some 
of  the  singly  diffracted  rays  which  they  produce.  One 
diffracted  ray  from  each  edge  is  shown  crossing  the 
slit  and  hitting  the  opposite  edge,  producing  doubly 
diffracted  rays. 
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produce  other  diffracted  rays.  These  rays  may  in  turn  produce  other  diffracted 
rays,  etc.  All  of  the  diffracted  rays — both  singly  and  multiply  diffracted  from 
edges  and  vertices — must  be  taken  into  account.  They  can  all  be  characterized 
by  a  more  general  extension  of  Fermat's  principle  which  is  similar  to  the  two 
special  forms  previously  described. 

One  or  more  diffracted  rays  pass  through  each  point  in  space.  These  rays 
account  for  the  field  in  the  shadow,  and  also  account  for  modifications  of  the 
field  outside  the  shadow.  As  we  shall  see, the  singly  diffracted  rays  are 
quantitatively  most  important  and  the  edge  diffracted  rays  more  so  than  the 
vertex  diffracted  ones. 

3.  The  field  on  a  ray 

Now  that  the  rays  have  been  characterized,  we  turn  to  the  quantitative 
description  of  the  field.  To  this  end  we  assume  that  there  is  an  amplitude  A 
and  a  phase  kT  associated  with  each  ray.  Both  of  these  quantities  may  vary  with  the 
distance  s  (from  some  origin)  along  the  ray.  In  terms  of  these  quantities  the 
field  on  the  ray  is  assumed  to  be  given  by  A(s)e      .  If  the  field  is  electro- 
magnetic, A  is  a  vector,  but  if  the  field  is  scalar—which  we  will  now  suppose — 
then  A  is  a  scalar.  In  neither  case  is  A  assumed  to  be  real. 

We  assume  that  the  difference  in  phase  between  two  points  on  the  same 
ray  is  equal  to  k  times  the  length  of  the  portion  of  the  ray  between  them.  Thus 
for  two  points  at  s  and  s  we  have  T(s)  -  T(s  )  ■  s  -  s  .  If  we  choose  s  ■  0 
and  denote  Y(0)  by  T  this  yields 


(2)         T(s)  -  TQ  ♦  s. 


We  also  assume  that  there  is  a  direction  of  propagation  on  each  ray  and  that  s 
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increases  in  this  direction.  The  direction  of  propagation  thus  determines  the 
direction  in  which  the  phase  increases.  From  these  considerations  we  see  that 
the  phase  can  be  found  at  any  point  s  on  a  ray  by  means  of  equation  (2)  provided  that 
it  is  known  at  some  other  point  on  the  same  ray  and  provided  that  the  direction  of 
propagation  is  alsc  known. 

To  determine  how  the  amplitude  varies  along  a  ray,  we  assume  that  the 
optical  form  of  the  principle  of  conservation  of  energy  is  satisfied.  This 
principle  states  that  the  flux  of  energy  is  the  same  at  every  cross  section  of  a 
tube  of  rays  in  a  non-dissipative  medium,  such  as  the  one  under  consideration. 

Furthermore  we  assume  that  the  flux  through  a  narrow  tube  of  cross  sectional  area 

2 

dcr  is  cA  do-  where  c  is  the  propagation  velocity.  By  applying  this  energy 

principle  to  two  cross  sections  of  such  a  tube,  we  obtain  the  relation 

2         2 

cA  do-  =  c  A  dcr  .  In  a  homogeneous  medium  the  propagation  velocity  is 

°  /^oY1-/2 

constant,  so  c  »  c  .  Thus  we  find  from  the  energy  relation  that  A  =  A  (  .■  ■  ■  ) 

If  we  imagine  that  the  tube  shrinks  down  upon  a  single  ray  then  this 

equation  relates  the  amplitudes  A  and  A  at  two  points  —  say  P  and  P  —  on  this 

ray.  In  this  relation  do-  /do-  is  the  limit  of  the  area  ratio  as  the  tube  shrinks 

to  the  ray.  This  limit  is  just  the  Jacob ian  of  the  mapping  of  one  cross  section 

of  the  tube  on  the  other  by  means  of  rays.  Thus  we  find  that  the  amplitude  A 

can  be  computed  at  any  point  P  on  a  ray  provided  the  amplitude  A  is  known  at 

some  other  point  P  on  the  same  ray.  To  perform  this  computation  it  is  not 

sufficient  to  deal  with  one  ray,  but  it  is  instead  necessary  to  take  account  of 

other  rays  near  the  one  under  consideration  in  order  to  determine  the  factor 

do-  /do-  . 
o 

The  area  ratio  dcr  /dcr  can  be  computed  very  easily  if  we  first  observe 
that  a  normal  cross  section  of  a  narrow  tube  of  rays  is  just  a  portion  of  a  wave- 
front.  Then  the  two  cross  sections  at  the  two  points  along  the  tube  are  portions 
of  parallel  wavefronts.  Now  the  ratio  of  areas  of  corresponding  small  portions 
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of  parallel  surfaces  is  the  reciprocal  of  the  ratio  of  the  Gaussian  curvatures 
of  these  surfaces.  To  compute  this  ratio  we  denote  by  p..  and  p?  the  principal 
radii  of  curvature  of  the  wavefront  at  the  point  P  on  the  ray.  Then  the  cor- 
responding radii  at  the  point  P  a  distance  s  along  the  ray  are  p,+s  and  p?+s. 
Thus,  since  the  Gaussian  curvature  is  the  reciprocal  of  the  product  of  the  prin- 
cipal  radii,  we  have 

(3)  d<ro  _  ,   P1P2 

K*}  do-   (P;l+sKp2+s;  ' 

Now  using  this  in  our  preceding  expression  for  A  we  obtain 

p1p2         \    1/2 


(U)  A  -  AQ 


[Pl+S]  [P2+S] 


In  order  to  combine  the  expressions   (2)  and  (U),  which  determine  the 
phase  and  amplitude  on  a  ray,  we  first  choose  P    as  the  origin  of  s.     Then  at 
the  point  s  on  the  ray  the  field  u  is  given  by 

($)  u  -    A, 


o 


plp2 


(p-j+sHpg+sj 


1/2e*<=*V 


From  (5)  we  see   that  for  large  s,  the  field  on  a  ray  decreases  like  s"  ( 

as  in  a  spherical  wave,  provided  that  both  p..  and  p?  are  finite.  If  one  radius 

-3/2 
is  infinite,  the  field  decays  like  s    , as  in  a  cylindrical  wave,  while  if  both 


*This  expression  for  do- /do- follows  at  once  from  the  geometrical  interpretation 
of  the  radii  of  curvature  shown  in  Figure  3.  Note  that  the  ratio  of  corres- 
ponding horizontal  lines  in  the  two  wavefronts  is  p^Ap^+s)  and  that  of  cor- 
responding vertical  lines  is  Pg/tPp+s).  The  ratio  of  the  areas  is  the  product 
of  these  two  ratios. 
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radii  are  infinite  the  field  is  of  constant  amplitude,  as  in  a  plane  wave.  In 
these  respective  cases  there  are  two,  one,  or  no  points  on  the  ray  at  which  the 
above  expression  for  the  field  becomes  infinite.   Such  points  are  called  caustic 
points  and  their  locus  is  the  caustic  associated  with  the  ray  system.  Since  the 
present  theory  yields  an  infinite  value  for  the  field  on  a  caustic,  it  is  obviously 
not  valid  at  such  a  surface.  Later  we  shall  see  how  this  defect  of  the  theory 
can  be  remedied. 

I,     Reflection  and  diffraction  coefficients 

The  preceding  considerations  have  resulted  in  equation  (5),  which  enables 
us  to  find  the  field  .&  any  point  on  a  ray  provided  that  it  is  known  at  one  point 
on  the  ray.  We  will  now  determine  the  field  at  the  point  of  initiation  of  each    ■ 
reflected  and  diffracted  ray,  assuming  that  it  is  known  on  the  corresponding 
incident  ray.  To  this  end  we  assume  that  the  initial  value  of  the  field  on  each 
reflected  or  diffracted  ray  is  proportional  to  the  value  of  the  field  on  the  corres- 
ponding incident  ray  at  the  point  of  reflection  or  diffraction.  The  respective 
coefficients  of  proportionality  will  be  called  reflection  and  diffraction  coefficients. 

Applying  these  considerations  to  the  field  u^  on  a  reflected  ray,  we 
have  at  the  point  of  reflection 


(6)        uref(0,y,z)  -  RHj^CO.y.a). 

In  (6)  R  is  the  reflection  coefficient.  This  coefficient  depends  upon  the  nature 
of  the  field  and  the  nature  of  the  screen  at  the  point  of  reflection.  It  may 
also  depend  upon  the  angle  of  incidence  and  upon  k.  As  an  example,  if  the  total 


.  vanishes  on  the  screen,  then  P  -  -1,  as  we  see  from  (6).  On 


xnc 


ref 
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the  other  hand,  if  the  normal  derivative  of  the  total  field  vanishes,  then  we 
find  that  R  =  +1.  Other  boundary  conditions  lead  to  other  values  of  R. 

Equation  (6)  determines  the  initial  value  of  the  field  on  each  reflected 
ray.  Therefore,  with  the  aid  of  equation  (5),  the  field  can  be  found  at  all 
points  on  each  reflected  ray.   To  calculate  this  field  it  is  necessary 
to  determine  the  radii  of  curvature  of  the  reflected  wavefronts  at  each  point 
of  reflection  on  the  screen.  Because  the  screen  is  plane,  these  radii  are  equal 
to  the  corresponding  radii  of  the  incident  wavefronts  at  each  point.  Therefore 
the  geometrical  factor  rpsPpApT+sHp^+s) I    and  the  phase  factor  e 
will  have  the  same  values  at  a  point  P  distant  s  from  the  screen  along  a  reflected 
ray  as  they  would  have  at  a  point  Q  the  same  distance  behind  the  screen  along 
the  corresponding  incident  ray   (see  Figure  U).  But  these  factors  would  give 
the  value  of  the  incident  field  at  Q.  Therefore,  since  the  y-  and  z-coordinates 
of  P  and  Q  are  the  same,  while  their  x-coordinates  differ  in  sign,  (5)  and  (6) 
yield 
^        uref^x,y,z^  "     R  uinc(-x,y,z). 

In  (7)  R  is  the  value  of  the  reflection  coefficient  appropriate  to  the  reflected 
ray  through  x,y,z.  If  no  reflected  ray  passes  through  x,y,z  then  u  «(x,y,z)  =  0. 

Let  us  now  consider  the  field  on  an  edge-diffracted  ray.  First  we  note 
that  neighboring  diffracted  rays  of  the  same  cone  intersect  each  other  at  the 
diffracting  edge.  Thus  the  edge  is  a  caustic  of  the  set  of  diffracted  rays 
and  each  point  of  it  is  a  center  of  curvature  of  the  diffracted  wavefronts. 
Therefore,  as  we  have  seen,  our  result  (£)  will  yield  an  infinite  value  for  the 
field  at  the  edge.  Consequently  we  cannot  determine  the  field  on  a  diffracted 
ray  simply  by  equating  the  diffracted  field  at  the  edge  to  a  multiple  of  the 
incident  field.  Instead  we  proceed  as  follows. 
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Consider  the  point  P    from  which  the  distance  s  is  measured  along  a 
diffracted  ray.     Let  p«  be  the  distance  from  PQ  to  the  edge  along  the  ray.     Then 
p.  is  one  radius  of  curvature  of  the  diffracted  wavefront  at  PQ,  because  one  center 
of  curvature  of  the  wavefront  is  at  the  edge.     Now  rewrite   (U)  in  the  form 


(8)  Ao  p1/2  -  A^sHp^s)]1/2  •    p-l/2 


If  the  right-hand  side  of  (8)   is  evaluated  at  some  fixed  point  P  on  the 
diffracted  ray,  then  it  is  independent  of  the  location  of  P  .     This  is  so 
becaus-3  A  is  a  constant  at  P  and  because  p,+s  and  p?+s  are  the  distances  from 
P  along  the  ray  to  the  two  centers  of  curvature  of  the  diffracted  wavefront. 
Therefore  the  left  side  of  (8)  is  also  a  constant  independent  of  the  location 

of  P  .     Let  us  denote  this  constant  by  A  .     Furthermore,  let  us  choose  P     at 

o  o  '  o 

the  edge.  Then  p„  ■  0  and  (8)  yields 

1/2 


(9)         A  "  A0 


rp^+sjs  j 


In  (9)  s  denotes  distance  from  the  edge  along  a  diffracted  ray. 

By  combining  (9)  and  (2)  we  obtain  for  the  field  u  .   diff  on  a  ray 
diffracted  from  an  edge,  the  result 


(10)        uedge  diff  "  Ao 


pl 


(ppsTs 


1/2  aiJc(s+T0) 


The  radius  of  curvature  p.  is  the  distance  from  the  edge  to  the  other  caustic  of 
the  diffracted  rays.  This  distance  is  computed  in  Appendix  I,  and  is  found  to 
be  given  by  (see  Figure  5) 
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Incident  Ray 


Figure  £ 
A   pair  of  neighboring  incident  rays  hitting  the 
curved  edge  of  a  thin  screen,  and  some  of  the 
diffracted  rays.  The  two  cones  of  diffracted 
rays  intersect  at  the  caustic,  which  is  at  the 
distance  p,  from  the  edge  along  one  of  the 
diffracted  rays. 


Diffracted 
Ray 


Figure  6 
The  projection  of  incident  and 
diffracted  rays  into  a  plane  normal 
to  the  edge.  The  angles  a  and  9  are 
the  angles  between  these  projections 
and  the  negative  x-axis,  measured  as 
shown  in  the  figure.  The  diffracting 
edge,  which  is  perpendicular  to  the 
plane  of  the  figure,  is  the  edge  of 
the  screen  or  the  edge  of  a  wedge 
(see  .Appendix  II). 
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(11)       p^^  --psin  p  (p  *p  sin  p  +  cos  6)   , 

In  (H),  P  <:  0  denotes  the  radius  of  curvature  of  the  edge,  p  is  the  angle 
between  the  incident  ray  and  the  (positive)  tangent  to  the  edge,  p  ia  the 
derivative  of  p  with  respect  to  arclength  along  the  edge, and  6  is  the  angle 
between  the  diffracted  ray  and  the  normal  to  the  edge.  It  is  to  be  noted  that 
equation  (ll)  is  independent  of  the  direction  in  which  arclength  increases 
along  the  edge  since  both  factors  p  and  sin  p  change  sign  when  this  direction 
is  reversed.  Further.nore  the  normal  to  the  edge,  which  lies  in  the  plane  of  the 
screen,  is  supposed  to  point  toward  the  center  of  curvature  of  the  edge. 

We  now  assume  that  the  factor  A  e   °,  which  determines  the  field  on  a 

o      * 

ikY 
diffracted  ray  according  to  (10),  is  proportional  to  the  incident  field  Ae 

at  the  point  of  diffraction.  We  will  denote  the  coefficient  of  proportionality 

by  D  and  call  it  the  diffraction  coefficient.  Thus  (10 )  becomes 

l/2 

<i2>       <w  «*  - D  A  [s<1^ls)l     eik(s*T>  • 

The  diffraction  coefficient  depends  upon  the  nature  of  the  field  and  the  nature 
of  the  screen  at  the  point  of  diffraction.  It  also  may  depend  upon  k  and  upon 

the  directions  of  the  incident  and  diffracted  rays.  It  is  to  be  noted  that 

-1/2 
u    diff  beh£Ves  liJce  a  cylindrical  wave  (s    )  near  the  edge,  and  like  a 

spherical  wave  (3"  )  far  from  the  edge,  provided  that  p.  is  finite. 

In  order  to  obtain  two  examples  of  diffraction  coefficients,  we  will  now 

consider  a  scalar  field  which  satisfies  the  reduced  wave  equation.  In  the  first 

case  we  suppose  that  the  field  vanishes  on  the  screen  and  in  the  second  case 

that  its  normal  derivative  vanishes  on  the  screen.  In  each  case  the  diffraction 

of  a  plane  wave  by  a  screen  in  the  form  of  a  half-plane  can  be  described  exactly 

by  the  solution  of  a  boundary-value  problem  for  the  reduced  wave  equation.  This 
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was   shown  by  Sommerfeld L l .  These  exact  solutions  can  also  be  expanded 
asymptotically  for  large  k.  On  the  other  hand,  these  problems  can  be 
treated  by  the  present  method,  using  (7)  for  the  reflected  field  and  (12)  for 
the  diffracted  field.  The  asymptotic  expansion  of  the  exact  solution  and  the 
field  computed  from  (7)  and  (12)  are  compared  in  Appendix  H  and  are  found  to 
agree  perfectly  provided  that 

in/U     r        -,  •,     i 

/,  ~s  D - sec  *(©-a)  +  esc  •*(»•♦« ). 

U3j  2  /m   sin  p  L   2  2    J 

The  angles  o  and  9  in  (13)  are  defined  as  follows:  Project  the 
incident  and  diffracted  rays  into  the  plane  normal  to  the  edge  at  the  point 
of  diffraction.  The  angles  between  these  projections  and  the  normal  to  the 
screen  are  a  and  6  respectively  (see  Figure  6 J  The  angle  p  is  the  angle  between 
the  incident  ray  and  the  edge.  The  upper  sign  in  (13)  applies  to  the  case  in 
which  the  field  vanishes  on  the  screen  and  the  lower  sign  to  the  case  in  which 
the  normal  derivative  vanishes  on  the  screen. 

Equation  (13)  gives  the  diffraction  coefficients  for  two  kinds  of  screens 
for  a  scalar  field  satisfying  the  reduced  wave  equation.  Thus  for  such  fields 
and  screens  the  field  on  an  edge  diffracted  ray  is,  from  (11),  (12)  and  (13), 

ik(Y+s)+ijiA  r   ,  ,     -, 

{lh)  uedee  diff  =  "  **!= lsec  2  (&><t)±  csc  ?  (e+<l) 

6(186  dU1    2  f£fc   sin  p   L    d  2     J 


;(1  m   3  [cos  £  +  pB  sin  p]  ^ 
p  sin  p 
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It  is  to  be  noted  that  D,  and  therefore  the  edge-diffracted  field  given 
by  (lii),  are  inversely  proportional  to  ^Tc   .     Thus  the  diffracted  field  decreases 
in  magnitude  as  k  increases,   i.e.   as  X  decreases.     This  is  in  accordance  with 


.  iu  - 


our  expectations  about  diffraction.  This  behavior  is  quite  general,  and  is  not 

restricted  to  the  special  fields  or  screens  which  we  have  just  considered.  It 

1/2 

follows  from  (12)  that  D  must  have  the  dimension  (length)  '  since  A  and 

u  ,   Hif-r  ^lave  "kne  same  dimensions.  The  only  length  upon  which  D  may  depend 

l/2       -l/2 

is  X  and  therefore  D  must  be  proportional  to  X    or  to  k    . 

Let  us  turn  now  to  the  corner-diffracted  rays.  Since  these  rays  emanate 
from  the  corner  in  all  directions,  the  corresponding  diffracted  wavefront3 
are  spheres  with  the  corner  as  center.  Both  centers  of  curvature  of  these 
wavefronts  are  at  the  corner,  which  is  consequently  a  caustic.  Considerations 
similar  to  those  which  led  from  (5)  to  (12)  in  the  case  of  edge  diffraction  now 
lead  from  (5)  to 

eik(s+T) 

(!£)         u™~^„  ji-r-e   =  C  A  = 
corner  diif      s 


ikY 
Here  s  denotes  distance  from  the  corner,  Ae     is  the  incident  field  at  the 

corner,  and  C  is  the  coefficient  of  corner  diffraction.         C,  like  D  and 

R,   depends  upon  the  nature  of  the  field,   the  nature  of  the  screen  at  the 

corner,  the  directions  of  the  incident  and  diffracted  rays,   and  the  value  of 

k.     In  addition  it  depends  upon  the  angle  between  the  two  edges  meeting  at  the 

corner.     From  (l5)  we  see  by  dimensional  considerations  that  C  must  be  proportionel 

to  X  or  k~  .     Thus  the  corner-diffracted  field  decreases  faster  than  the  edge - 

diffracted  field  as  k  increases. 

In  order  to  determine  C  for  a  scalar  field  satisfying  the  reduced  wave 

equation  and  for  screens  on  which  the  field  or  its  normal  derivative  vanishes, 

we  may  proceed  as  follows.     Consider  a  plane  wave  incident  upon  an  infinite 

corner,   i.e.   a  screen  with  two  semi-infinite  straight  edges  meeting  at  a 

corner.     Solve  this  boundary-value  problem  exactly  for  the  reduced  wave 

equation  and  expand  the  solution  asymptotically  for  large  k.       Compare 
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this  result  with  the  field  given  by  (7),  (lU)  and  (l5).  The  two  results  should 
agree  for  an  appropriate  choice  of  C,  which  will  then  be  the  desired  corner- 
diffraction  coefficient.  This  procedure  has  not  yet  been  carried  out  because 
of  the  intractable  form  of  the  solution  of  the  corner-diffraction  problem,  which 
has  been  considered  by  L.  Kraus  "■  ■• . 

5.  The  total  field 

We  have  now  determined  all  the  rays,  found  how  the  field  varies  along 
each  ray  and  determined  tha  'initial1  value  of  the  field  on  each  reflected  or 
singly  diffracted  ray.  Exactly  the  same  formulas  apply  to  a  doubly  diffracted 
ray,  the  corresponding  singly  diffracted  ray  being  considered  as  W  incident 
ray.  In  the  same  way  the  field  on  every  multiply  diffracted  ray  can  be  deJ  rmined. 
We  are  now  ready  to  consider  the  total  field. 

We  assume  that  the  total  field  u  at  any  point  P  is  the  sum  of  the  fields 

on  all  the  rays  through  P.  We  also  define  the  diffracted  field  at  P,  u, .««,  as 

the  sum  of  the  fields  on  all  diffracted  rays  through  P.  Finally  we  define  the 

geometrical  optics  field  at  P,  u    =  u.   +  u  „  as  the  sum  of  the  fields  on 

geom    inc    ref 

all  the  incident  and  reflected  rays  through  P.  Then  we  can  write 


(16)        u  =  u    +  u, .„„ 

geom    diff 

The  geometrical  optics  field  is  zero  in  the  shadow  region,  through  which 

neither  incident  nor  reflected  rays  pass.  Thus  in  the  shadow  the  total  field 

is  "U^f**  However  u,-»„  is  also  present  in  the  'illuminated'  region,  in  which 

u  ^  is  not  zero, 
geom 

The  diffracted  field  u , . „„  consists  of  a  finite  or  infinite  sum  of  the  fields 

ail  I 

on  all  the  singly,  doubly  and  multiply  diffracted  rays  through  each  point.  These 
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rays  may  be  diffracted  from  edges,  from  corners,  or  from  both.  Each  time  a  ray 

is  diffracted  by  an  edge,  its  field  is  multiplied  by  a  diffraction  coefficient  D 

-1/2 
which  contains  the  factor  k    .  Each  time  it  is  diffracted  by  a  corner,  its 

field  is  multiplied  by  a  corner  diffraction  coefficient  C  which  contains  the 

factor  k~   .     Thus  for  large  values  of  k,   the  fields  on  multiply  diffracted  rays 

become  rapidly  weaker  as  the  number  of  diffractions  increases.  Pbr  very  large 
values  of  k,   the  singly  edge-diffracted  field,  u  in_-i,,  eA„e  diff'   ^s  mos-t 

significant,   all  the  other  diffracted  fields  being  much  weaker.  Even  the 

singly  corner-diffracted  field  is  insignificant  compared  to  it.  Of  course, 

in  those  regions  of  the  shadow,   if  any,  where  ugine-ie  eH2e  diff  ^3  zer0» 

u   .     .  . „„  becomes  important, 

single  corner  diff 

The  singly  edge -diffracted  field  at  any  point  P  can  be  written  as 

<«>  Vngla  edge  di^»"    Z  B  A  [.(l^lfV''*", 

The  sum  extends  over  all  the  singly  edge-diffracted  rays  through  P.     In  case  D 
is  given  by  (13 )y   (17)  becomes 

.   ik(Y+s)+inA  r         ,  ,1 


2  yzitk"  sin  p 


-1/2 


Cl        s[cos  6  +  pp  sin  p]  \ 
p  sin  p 


The  formula  (18)  expresses  the  singly  diffracted  field  explicitly  in  terms  of 
geometrical  quantities.     This  formula  will  be  applied  in  the  following  sections 
to  apertures  and  disks  of  various  special  shapes.     In  certain  cases  the  corres- 
ponding formula  which  includes  all  multiply  diffracted  rays  will  also  be  used. 
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Before  we  apply  our  result  to  special  cases,  we  will  make  certain 
general  comments  concerning  it.  In  the  first  place,  it  is  clear  from  the 
derivation  that  the  result  (12)  applies  to  diffraction  by  any  edge,  and  not 
merely  by  the  edge  of  a  thin  screen.  For  example,  it  applies  to  diffraction  by 
the  edge  of  a  wedge,  provided  that  the  appropriate  diffraction  coefficient  is 
used.   Ihis  coefficient  is  deduced  for  two  kinds  of  wedges  —  those  on  which  the 
field  vanishes  and  those  on  which  its  normal  derivative  vanishes  —  in  Appendix  II. 
Similarly,  equation  (l£)  applies  to  diffraction  by  any  corner  or  vertex,  and  not 
merely  to  a  corner  on  the  edge  of  a  thin  screen.  Thus  it  applies  to  diffraction 
by  the  tip  of  a  cone  or  the  corner  of  a  cube,  provided  that  the  appropriate 
diffraction  coefficient  is  used  in  each  case. 

Next  it  is  important  to  observe  that  the  diffraction  coefficients  can 
always  be  obtained  from  certain  exact  special  or  canonical  solutions  of  boundary- 
value  problems  for  the  differential  equations  governing  the  field.  The  geometry  in 
the  canonical  problem  may  be  very  simple,  provided  only  that  it  coincides  locally  with 
the  geometry  near  the  point  of  diffraction.  Therefore  it  may  involve  just  a 
plane  incident  wave  -and  a  wedge-like  or  conical  obstacle. 

Canonical  solutions  can  also  be  used  to  correct  the  results  of  the 
present  theory  at  and  near  caustics  and  shadow  boundaries,  where  the  present  theory 
yields  infinite  values  for  the  field.  Thus,  for  example,  at  an  edge  the  diffracted 
field  is  infinite.  Furthermore  the  geometrical  optics  field  is  discontinuous 
across  the  shadow  boundary  which  passes  through  the  edge.  Near  the  edge  these 
two  fields  may  be  replaced  by  the  field  obtained  from  the  solution  of  a  boundary- 
value  problem  for  diffraction  by  a  half-plane.  The  edge  of  the  half -plane  must  be 
tangent  to  the  actual  aperture  edge.  Of  course  the  incident  field  and  the 
properties  of  the  screen  in  the  half-plane  problem  must  correspond  to  those  in 
the  actual  problem.  This  procedure  is  just  the  one  which  was  proposed  by  Braunbek^-" 
as  a  basis  for  the  calculation  of  diffraction  by  an  aperture  in  a  thin  screen. 
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In  Part  IT  of  this  article  we  will  show  that  his  proposal  can  be  worked 
out  in  detail  for  an  aperture  of  any  shape  and  that  it  leads  exactly  to  our 
equation  (18).  It  then  follows  that  our  method  is  equivalent  to  Braunbek's  as 
far  as  rays  singly  diffracted  from  edges  are  concerned.  However  his  methods 
does  not  t^ke  account  of  multiple  diffraction,  nor  does  it  yield  the  correct 
field  on  corner-difu.x'acted  rays.  Furthermore,  it  requires  an  integration  in 
order  to  compute  the  field,  while  the  present  method  does  not. 

In  Part  II  we  will  also  consider  the  various  forms  of  the  Kirchhoff 
method  and  show  that  each  of  them  leads  to  a  result  of  the  form  given  by 
equation  (17).  However  the  diffraction  coefficients  D  will  not  coincide  with 
the  exact  diffraction  coefficients  given  by  equation  (13) • 

It  should  be  observed  that  the  diffraction  coefficients  given  by 
equation  (13)  become  infinite  on  the  shadow  boundary  ©  =  -a  and  on  the  boundary 
of  the  reflected  field  6  =  n   +  a.  These  are  the  surfaces  across  which  the 
geometrical  optics  field  u.   +  u  »  is  discontinuous.  Thus  the  present  results 
are  invalid  on  these  surfaces. 

The  modifications  of  the  preceding  results  which  must  be  made  in  the 
case  of  electromagnetic  fields  are  rather  minor.  If  the  field  u  denotes  an 
electric  or  a  magnetic  field  then  it  must  be  a  vector.  We  assume  that  all  of 
our  preceding  results  apply  to  these  fields  when  we  consider  the  amplitude  A  to 
be  a  vector  whose  direction  remains  constant  along  a  ray  and  which  is  perpendicular 
to  it.  If  we  consider  both  electric  and  magnetic  fields,  they  must  also  be 
perpendicular  to  each  other.  Since  A  io  a  vector,  the  reflection  coefficient  R 
and  the  diffraction  coefficients  D  and  C  must  now  be  matrices.  The  matrices  R 
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for  reflection  of  the  electric  or  magnetic  field  from  a  perfect  conductor 
or  from  a  dielectric  are  known.  These  matrices  are  given  by  the  well-known 
exact  solutions  of  Maxwell's  equations  for  the  problems  of  reflection  of  a 
plane  electromagnetic  wave  from  a  plane  conductor  and  from  a  plane  dielectric 
interface.  Consequently  they  will  not  be  given  here.  However  in  Appendix  III 
we  consider  the  diffraction  of  an  electromagnetic  wave  by  the  edge  of  a  per- 
fectly conducting  screen  and  determine  the  diffraction  matrix  D. 

We  will  now  apply  our  results  to  various  problems,  at  first  assuming 
the  field  to  be  scalar  and  the  diffraction  coefficient  D  to  be  given  by  (13). 


6.  Diffraction  of  spherical  and  cylindrical  waves  by  a  half-plane 

As  a  first  application  of  our  result,  let  us  find  the  field  produced 

when  an  arbitrary  wave  is  incident  upon  a  half-plane.  In  this  case  no  multiply 

diffracted  or  corner-diffracted  rays  occur.  We  will  assume  that  the  angle 

between  the  incident  ray  and  the  edge  varies  monotonically  along  the  edge, 

as  is  the  case  with  incident  plane,  spherical  and  cylindrical  waves.  Then 

there  is  just  one  diffracted  ray  through  each  point  not  on  the  edge.  Thus 

u, .ff  consists  of  a  single  term  which  is  given  by  equation  (12)  or  (lit).  These 

formulas  can  be  simplified  since  p,  the  radius  of  curvature  of  the  edge,  is 

infinite.  It  then  follows  from  (ll)  that  p-  =  -  p   sin  p. 

If  the  incident  wave  is  cylindrical  and  if  its  axis  is  parallel  to  the 

edge,  then  (3  =  n/2  for  all  rays  incident  on  the  edge.  The  field  in  a  cylindrical 

ikr 
wave  is  of  the  form  a ,  where  r  denotes  distance  from  the  axis  and  a^  is  a 

0  y&  ° 
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constant.  This  follows  from  arguments  similar  to  those  used  in  deriving 
equation  (12).  We  choose  a  =  /—  e~    ,  so  that  the  cylindrical  wave  is  that 
due  to  a  line  source  of  unit  strength.  Now  if  r  is  the  distance  of  the  line 
source  from  the  edge,  then  upon  inserting  the  above  expression  for  the  cylindrical 
wave  into  (lU)  we  obtain 


(19) 


,ik(s+r0) 


u 


diff 


27ik)^s" 


sec  i  (e-a)  ±  esc  -|  ($+a)  , 


Equation  (l9)gives  the  field  due  to  diffraction  of  a  cylindrical  wave  from  a  screen 
on  which  u  vanishes,  if  the  upper  sign  is  taken,  or  on  which  the  normal  derivative 
of  u  vanishes,   if  the  lower  sign  is  taken.     The  boundary-value  problem  for  the 
reduced  wave  equation  can  be  solved  exactly  for  such  screens  with  a  cylindrical 
incident  wave,   and  the     solutions  can  be  expanded  asymptotically  for  large  k. 
The  results  of  these  expansions  are  precisely  given  by  equation  (19)  plus  u 
This  agreement  constitutes  a  check  on  our  theory. 

Next  let  us  suppose  that  the  incident  wave  is  spherical,  with  its  center 

at  a  point  0.     Then  if  r  denotes  distance  from  0,  the  incident  field  is  of  the 

ikr  / 
form  a  e       It  and  we  will  set  a     -  1.     This  form  of  the  field  can  be  obtained 
o  o 

by  arguments  like  those  which  lead  to  (1$).     At  a  point  Q  on  the  edge  at  distance 

•-1       / 
r  from  0,  p  --r/sinp  where  p  is  the  angle  between  the  ray  from  0  to  Q.     Since 

p,   "-P™  sin  p,  we  see  that  p,   *  r  at  Q.     Using  these  results  in  (lU)  we  obtain 


(20) 


ik(r+s)+inA 


u 


diff 


2  sin  p  /2~nk  /rs(r+s) 


sec  -x  (6-a)  +  esc 


\  (*•>]  . 


Equation  (20)  gives  the  field  due  to  diffraction  of  a  spherical  wave  from  a  screen 
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on  which  u  or  its  normal  derivative  vanishes,  according  as  the  upper  or  lower 

sign  is  taken.     The  boundary-value  problem  for  the  reduced  wave  equation  has 

been  solved  exactly  for  such  screens  with  a  spherical  wave  incident.     The 

solution  can    be     expanded  asymptotically  for  large  k  and  the  result  agrees 

precisely  with  u, .„„  in  equation  (20)  plus  u         .     This  agreement  is  another 
ail 1  geom 

check  on  our  theory. 

7.     Diffraction  by  a  slit 

As  a  second  application  of  our  result  we  will  compute  the  field  produced 
by  a  plane  wave  incident  upon  a  slit  of  width  2a  and  of  infinite  length,   in  a 
thin  screen  (see  Figure  2).     We  will  suppose  that  each  incident  ray  lies  in  a 
plane  normal  to  the  edges  of  the  slit,  so  that  3  =  n/2  for  all  rays  incident 
upon  the  edges.     This  problem  is  essentially  two-dimensional,   since  each 
diffracted  ray  will  also  lie  in  a  plane  normal  to  the  edges.     The  more  general 
case  in  which  p  4  n/2  can  be  treated  almost  as  easily,  or  the  results  for  it  can 
be  obtained  from  the  present  ones  by  the  method  explained  in  Appendix  II. 

As  we  have  seen,   two  diffracted  rays  of  each  order  will  pass  through 
each  point.     To  compute  the  field  on  the  diffracted  rays  we  make  use  of  the  facts 
that  the  radius  of  curvature  of  each  edge  is  infinite  and  that  P  -  n/2, so  P  ■  0. 
It  then  follows  that  p..    is  also  infinite.     Now  using  (18)  we  can  write  the  singly 
diffracted  field  provided  that  we  know  the  incident  field.     Let  us  assume  that 
the  incident  field  is  e^x  cos  a  "  y  sin  a).     Then  (18)  yields 

ik(r1-a  sin  a)+in/l|.   r         ..  ,  -j 

(21)        "single  diff  "  -    2yCTF- [SeC  *  (<*1+  G)  -  CSC  *  (91  "  0)J 


ik(r2  +  a  sin  a)  +  in/fy. 


2  /2nkrr 


sec  ^  (62-  a)  ±  esc  ^  (62+  a)  . 
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In  equation  (21)  r,  and  r_  denote  the  distances  of  a  point  from  the  upper  edge 
at  (0,  +a)  and  from  the  lower  edge  at  (0,  -a)  respectively.  The  angles  &,  and 
9p  are  determined  by  the  corresponding  rays,  as  shown  in  Figure  2. 

Before  computing  the  fields  on  the  multiply  diffracted  rays,  we  will  exam- 
ine the  singly  diffracted  field  far  from  the  slit.  To  this  end  we  consider  the 
polar  coordinates  r,  0  of  a  point  in  the  x,y-plane.  If  this  point  is  fr.r  from 
the  slit  then  we  see  that  r,  ~  r  -  a  sin  0,  r_  ~  r  +  a  sin  0,  ©,  *>  n   +  0,  and 
6«  ^  n   ~  0*     Upon  using  these  expressions  in  equation  (21)  we  obtain 


<«>        Si^udiff-  -&>**"**  *.»> 


In  equation  (22)  we  have  introduced  f  ($),   the  radiation  or  diffraction  pattern 

s 

of  the  slit  due  to  singly  diffracted  rays.     From  (21)   it  is  found  to  be   given  by 

(??)  f  (d)  3  ■?   sin   [ka  (sin  0  *  sin  a)j        cos   [ka  (sin  0  +  sin  a)j 

s  k  sin  |  (0  +  a)  k  cos  I  (0  -  a) 


In  Figures  7  and  8,    |k  f  (ff)|    is  shown  as  a  function  of  0  for  normal 
incidence  (a  ■  0).     The  angle  0  ranges  between  0  and  n/2,   corresponding  to  the 
region  behind  the  slit.     Figure  7  is  based  upon  the  value  ka  *  8,   corresponding  to  a 
slit  about  3  wavelengths  wide,   and  Figure  8  corresponds  to  a  slit  10  wavele:igth» 

wide,   for  which  ka  ■  lOn.     It  is  to  be  noted  that  the  number  of  maxima  in  each 
curve  equals   the  number  of  wavelengths  which  fit  into  the  slit.      Furthermore  it 
is  interesting  that  the  curves  are  the  same  for  both  the  upper  and  lower  signs  in 
(23),   i.e.   for  both  the  screen  on  which  u  =  0   and  that  on  which  the  normal 
derivative  of  u     is  zero. 

For  comparison  the   curve  based  upon  the  exact  solution  of  the  boundary- 
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value  problem  for  the  reduced  wave  equation  is  also  shown  in  Figure  7.  This 

curve  applies  to  the  boundary  condition  u  «  0  on  the  screen. 

According  to  the  cross  section  theorem,   the  transmission  cross  section 

of  the  slit  (per  unit  length)  is  equal  to  tae  imaginary  part  of  the  diffraction 

pattern  in  the  forward  direction  0  =  -  a.  Applying  this  theorem  to  f  (0)  we 

s 

find  that  the  transmission  cross  section  is  2a  cos  a.  Since  this  is  the 

projection  of  the  slit  width  on  the  direction  normal  to  the  incident  rays,  it 

is  just  the  geometrical  optics  cross  section.  Thus  the  singly  diffracted  rays 

do  not  contribute  to  the  cross  section  when  it  is  computed  from  f  (0)  by  the 

cross  section  theorem.  Of  course  they  would  contribute  if  the  cross  section  were 

computed  by  integration  of  f  (0)  with  respect  to  0.  This  discrepancy  does  not 

s 

contradict  the  cross  section  theorem  because  the  theorem  applies  only  to  the 

exact  solution  of  the  boundary-value  problem  for  the  reduced  wave  equation. 

The  above  calculation  of  the  cross  section  did  not  take  explicit  account 

of  u nor  of  the  fact  that  u  .  ,  J.*..*,  becomes  infinite  on  the  two  boundaries 

geom  single  diff 

of  the  shadow.  Actually  u    ,  which  is  discontinuous  on  the  shadow  boundaries, 

geom'  ■ 

should  merge  with  us.j„p.-ie  ^iff  as  r  increases.  This  merging  is  automatically 

accounted  for  in  the  calculation  of  u  .  ,  ,.~„   for  r  large.  In  this  calcula- 

single  nl.ii 

tion  the  singularities  on  the  two  shadow  boundaries  cancel  each  other  and  the 
resulting  finite  expression  includes  the  effect  of  u    .  This  is  shown  by  the 
cross  section  result.  However  in  other  problems  the  same  simple  calculations  do 
not  suffice  and  other  considerations  are  needed.  This  will  be  seen  in  the  case  of 
the  circular  aperture. 

In  order  to  obtain  a  first  correction  to  the  geometrical  optics  cross 
section  we  need  merely  take  account  of  the  two  doubly  diffracted  rays  which  go  in 
the  direction  0  =  -  a.  To  do  this  we  note  that  the  singly  diffracted  ray  which 
leaves  the  lower  edge  in  the  direction  ©«  ■  n/2  hits  the  upper  edge.  It  reaches 


the 


upper  edge  when  rp  ■  2a.  Thus  the  field  on  this  ray  at  the  upper  edge  is, 
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if  we  choose  the  upper  sign  in  (21), 

ik(2a  +  sin  a)  +  in  A     , 
(21*)      -  2 sec  ±  (£  -  e) 

2  /SE5  2  2 

If  we  had  chosen  the  lower  sign  in  (21)  rather  than  the  upper  one,  we  would 
have  found  that  the  field  on  the  ray  &„  ~   n/2  is  zero.  Therefore  the  field  on 
all  doubly  and  multiply  diffracted  rays  would  be  zero  in  this  case,  and  (22) 
would  give  the  complete  diffracted  field.  We  will  reconsider  this  case  later, 
but  now  we  will  proceed  with  our  calculation  for  the  screen  on  which  u  «  0. 

To  find  the  field  on  the  doubly  diffracted  rays  produced  by  diffraction 

at  the  upper  edge,  we  again  use  equation  (lU),  treating  (2U)  as  the  incident  field. 

We  proceed  in  a  similar  way  to  find  the  field  on  the  doubly  diffracted  rays 

produced  at  the  lower  edge.  Combining  the  fields  on  the  two  doubly  diffracted 

rays  in  the  direction  0   ■  -  a  we  obtain 

(25)     u(r,-a)     +  tt(r,-o)      -  -  JS-  •lkpfi,lAffJI(-e)  ♦  f.(-a)}  . 
single  diff    double  diff    '   r         L        a 

Here  f  .(-a)  is  defined  by 
a 


26)     f,(-a) 


eiJ2ka(l+sin  a)  +  nA]       ^i[2ka(l-sin  a)  +  ttA] 


d  k/nH   I  1  -  sin  a 


We  now  apply  the  cross  section  theorem,  which  requires  that  we  take  the 
imaginary  part  of  f  (-a)  +  iJ-c).     We  thus  obtain  for  the  transmission  cross 
section  o~    the  expression 

,_„»                        _                            1        f  cos    [2ka(l+sin  a)  -  nAl        cos    [2ka(l-sin  a)   -  n/h] 
(27)  c    -  2a  cos  a  -    ■ ■ t— *•  +  ■ ' — l 


k/nka 


1  +  sin  a  1  -  3in  a 

-3/2 
The  correction  to  the  geometrical  optics  cross  section  is  of  the  order  k  *   .  This 
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result  applies  to  a  slit  in  a  screen  on  which  u  ■  0.  The  result  (2?)  has  also  been 

obtained  by  H.  Levinet9l#  jn   -the  special  case  of  normal  incidence, a  ■  °  and  (27) 


becomes 


(28)      ersl_cos(2ka-*/>4) 
2a       fH   (ka)3/2^ 


A  graph  of  cr/2a  as  a  function  of  ka  is  shown  in  Figure  9,   based  upon 
equation  (28).  For  comparison  a  curve  based  upon  the  exact  solution  of  the 
boundary-value  problem  for  the  reduced  wave  equation  is  also  shown. 

We  will  now  determine  the  field  on  each  multiply  diffracted  ray.  Then  we 
will  add  all  these  fields  together  to  obtain  the  total  multiply  diffracted  field. 
For  the  reason  previously  stated,  only  the  case  in  which  u  ■  0  on  the  screen  will 
be  considered* 

The  different  kinds  of  multiply  diffracted  rays  are  classified  in  the  fol- 
lowing table : 


Type 

Initial 
Diffraction 

Number  of 

at  each 
Lower 

diffractions 
edge 

Upper 

Final 
Diffraction 

A 

Lower 

n  +  1 

n 

Lower 

B 

Lower 

n 

n 

Upper 

C 

Upper 

n 

n  +  1 

Upper 

D 

Upper 

n 

n 

Lower 

The  first  column  designates  the  type  of  ray.  The  second  column  states  which  edge  is 
first  hit  by  the  incident  ray  that  gives  rise  to  a  particular  multiply  diffracted  ray. 
The  last  column  states  the  edge  at  which  the  multiply  diffracted  ray  is  finally  dif- 
fracted. The  other  columns  give  the  numbers  of  times  the  ray  is  diffracted  at  each 
edge.  The  integer  n  is  at  least  unit%  since  the  rays  are  all  multiply  diffracted. 

To  determine  the  field  on  any  ray  we  make  repeated  use  of  equation  (±h)» 
The  calculation  is  simplified  if  we  introduce  the  function  F(s,a,8)  defined  by 
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iks+in/U        r         ,  ,  -r 

(29)  F(s,a,e)     -     -  sec  i  (©-a)  ♦  csc  4  (9+a)       • 

2  /2lkT  L        '  £  J 


Then  according  to  equation  (1U),  the  diffracted  field  in  the  present  problem  is 
simply  the  incident  field  at  the  edge  multiplied  by  F(s,a,©).     Thus,  for  example, 
the  field  (2k)  at  the  upper  edge  on  the  ray  singly  diffracted  at  the  lower  edge  is 
just  e  a  P(2a,a,j.)c     By  referring  to  the  table  above,  one  easily  finds  that 

the   fields  on  the  various  rays  are  given  by 

A)  •**  "in  °  F(2.,a,J)  F2n"1(2a,J,  J)  F<r2£,  ©2), 

B)  •*•  8in  a  F(2a,a,J)  F2n-2(2.,J  ,  £)  FCr^J  ,  0^, 
C;  e"ika  ,in  <XF(2a,-a,5>  F2n-1(2.,J,  JWr^,©^, 
D)  eika  8in  »F(2at-«,J)  F20-2^,^)*^,  ©2). 

The  total  Multiply  diffracted  field  u^^  dlff  *»  now  obtained  by  adding 


(30) 


together  the  fields  on  the  four  types  of  rays  and  summing  with  respect  to  n. 

c 
The  only  series  to  be  summed  is  the  geometric  series 

r      2  "5-1  * 

which  has  the  sum  jl-F  (2a,*/2,H/2)j   .  Thus  we  find 

(31)     umult  diff  "  \   ,lka  8lD  tt*(2«.a,«/2)  [F(2a,«/2,«/2)  F  (r2,n/2,©2)  ♦  F^nA,©^] 

<¥(2a,-a,m/2)[F(2a,n/2,ii/2)F(r1,ii/2,01)+F(r2,ii/2,©2)j 


-ika  sin 


-1 


1  -  F2(2a,«/2,«/2) 


The  complete  diffracted  field  u^  -  ugingle  d±ff  +  u^  diff  is  given  by 
equations  (21)  and  (31)*  Exactly  the  same  result  has  been  obtained  by  S.N»  Karp 
and  A.  Russek"  '  by  a  very  different  method. 


-  27  - 


The  expression  (31)  for  the  multiply  diffracted  field  can  be  consid- 
erably simplified  at  points  far  from  the  slit.  By  using  the  notation  and 
approximations  explained  just  preceding  equation  (22),  we  find  that  (31) 
simplifies  to 

(32)     u  „  „„,  ■  -  /K  °lkr*i',/U  *L<«. 


m 


%ult  diff  "  "  /2nr  '       "m' 
The  radiation  or  diffraction  pattern  f„(0)  due  to  multiple  diffraction  is  given  by 

ika(sin  a  +  sin  0) 


(33)     fm(?0  »  [UnUcae-^M]-1   -  |  [sec  (f  -  |)  sec  (|  -  f ) 


e 


/it   o»     /it  .  0\  -ika(sin  a  +  sin  0)1 
+  sec  (j-  +  ^)  sec  (£  +  *)e  r  I 

.  0  /na  -i2ka-i£  [  M  ,n      o,     ,n   0s  ika(sin  a  -  sin  0) 
+  2  /  -s—  e      •*    sec  vj-  -  -^J  sec  (.-j-  +  £;e 


,n   t  a%     /n   0\  -ika(sin  a  -  sin  0) 
+  sec  Or  +  -j)  sec  (£  -  £)e 


when  a  ■  0  equation  (33)  simplifies  to 
(3U) 


fm(0)  -  -  k   r  yCnES  e^2*3"1*  +  l]"1  cos  §  C03(ka  si"  0)-i  sin  g  sin(kasin0) 

cos  0 

A  graph  of  k  times  the  total  diffraction  pattern,  k  <f  (0)  +  f  (0)}   ,  is 
shown  in  Figure  7  for  normal  incidence  (a  -  0)  with  ka  -  3*. 

To  find  the  transmission  cross  section  of  the  slit,  we  again  apply  the 
cross  section  theorem.  Now  we  must  compute  the  imaginary  part  of  fg(-«)  +  ;fm(~<t)• 
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In  the  case  of  normal  incidence  (a  ■  0)  this  yields 

_      cos  (2ka  -  y )  r        i/o  -,1-1 

(3$)        ~  -  1  -  — 575-2-  [i  -  (nka)-1/2sin(2ka  -  jHCUnka)"1! 

A  graph  of  o~/2a  as  a  function  of  ka  for  normal  incidence  (a=0)  is  shown  in 

Figure  9>  based  upon  equation  (35).  Corresponding  curves  based  upon  equation  (28) 

and  upon  an  exact  solution  of  the  reduced  wave  equation  are  also  shown.  It  can 

be  seen  that  both  the  result  (28)  based  upon  double  diffraction  and  equation  (3!>) 

are  quite  close  to  the  exact  curve.  However,  the  curve  based  upon  (35)  is  closer, 

especially  for  ka  small,  where  multiple  diffraction  is  important. 

From  (33)  we  see  that  f  (00  contains  as  a  factor  the  reciprocal  of 

m 

(Unika  e"    +  l).  If  this  quantity  is  zero,  then  f  (#)  and  the  field  u  will 
not  be  defined.  The  values  of  k  for  which  this  occurs  are  all  complex,  and 
correspond  to  the  free  vibrations  or  modes  of  the  field.  These  free  vibrations 
or  modes  contain  no  incident  waves,  but  consist  entirely  of  outgoing  waves.  Thus 
they  must  decay  exponentially  in  time,  and  we  see  that  this  is  so  from  their 
angular  frequencies  co  =  kc.  The  corresponding  factor  in  (3U)  is  the  reciprocal 
of  (ylmka  e"    "  *  +  1).  This  quantity  is  a  factor  of  the  previous  one,  and  thus 
contains  only  gome  of  the  roots.  These  roots  must  correspond  to  modes  symmetric 
about  the  x-axis  since  they  occur  in  the  case  of  normal  incidence.  The  other 
modes  are  consequently  antisymmetric.  All  the  roots  of  this  equation  have  negative 
imaginary  parts,  corresponding  to  decaying  modes. 

8.  Diffraction  by  a  circular  aperture 

Wis  will  now  apply  our  results  to  diffraction  of  a  normally  incident 
plane  wave  by  a  circular  aperture.  In  this  case  two  singly  diffracted  rays  pass 
through  each  point  which  is  neither  on  the  edge  nor  on  the  axis  of  the  aperture. 
They  come  to  any  point  P  from  the  nearest  and  farthest  points  on  the  edge.  Thus 
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the  singly  diffracted  field  given  by  (18)  contains  just  two  terms.  To  evaluate 
them  we  note  that  0  ■  n/2  for  all  incident  rays,  so  0  -  0.  If  the  radius  of 
the  hole  is  a,  then  p  ■  a.  Assuming  that  the  incident  wave  is  e3"^  we  obtain 
from  (18) 


eikr1+inA  [    &i      91  1 T      1 
single  diff(P)  "  "  T-= |_sec  T  *  csc  T      rl(l-a"  rlsin  91) 


u_ 

2  /2"r3c 

(36) 


-1-1/2 


ikr2+inA  r       eo  »,ir  i  1-1/2 


2  ys& 


r    2       2  1 P     -1 

lsec  -j-  +  csc  y   r2 (l-a  r2sin  ©g) 


In  (36)  r,  and  r«  denote  the  distances  from  P  to  the  two  points  of  diffraction 
on  the  edge,  while  9,  and  ©„  are  the  angles  determined  by  the  corresponding  rays, 
In  writing  (36)  we  used  the  relations  6,  =9.  -  ■*  and  6„  ■  9?  -  £. 

To  simplify  (36)  we  consider  a  point  P  far  from  the  hole  and  introduce 
the  polar  coordinates  r,  0  in  the  plane  containing  P  and  the  axis.  In  this 
plane  the  ray  configuration  and  the  angles  are  the  same  as  in  the  case  of 
normal  incidence  on  a  slit,  which  is  shown  in  Figure  2.  Therefore  the  approxi- 
mations made  in  that  case,  preceding  equation  (22),  are  applicable  again  and 
(36)  becomes 


ikr 
(37)     Vngla  dlff  "-  ^~  *.<*>• 

The  radiation  or  diffraction  pattern  f  (0)  due  to  single  diffraction  is  given  by 

s 


-  30  - 


(38) 


f.<W 


sin(ka  sin  0   -  |)    cos(ka  sin  0  -  •£) 


sin 


cos 


A  graph  of  |kf  (0)|   as  a  function  of  0,  based  upon  (38),  is  shown  in 
s 

Figure  10  for  ka  -  3n,     Ps  we  see  from  (38)  and  from  the  graph,  f  (0)  becomes 
infinite  at  0  =  0.     This  is  because  the  axis  0  «  0  is  a  caustic  of  the  singly- 
diffracted  rays  and  because  the  shadow  boundary  r  =  a  also  extends  in  the  direc- 
tion 0  »  0. 

Let  us  now  examine  the  field  (36)  at  points  near  the  axis.     We  set 

?     2l/2 
p  =  (y  +z  )  '     and  consider  points  for  which  p  is  small.     For  such  points 

we  set  ^  =  I  +  61  and  9g  =  |  +  62  where  6±  -  tanSc/U-p)  and  62  =  tan^-x/U+p). 

Then  r,   ■  x  sin  6.  +  (a  -  p)cos  5.  and  r«  =  x  sin  6_+  (a+p)cos  6*.    Thus 

1  -  a"  r,sin  9,  =    a"  p  and  1  -  a"  r2sin  &2    ■   -a"  p.     Finally  we  note  that 

rl~  r2*'(x2+a2)1^2  and  that  6i~  62^6"  tan"1x/a.     Inserting  these  results  into  (36)yieMs 


(39) 


usingle  diff 


v-     ik(x  sin  6  +  a  cos  6)    r         _  c  1 

/a  e  /6   ,   tin  ,0  .   ii\ 


•   cos(kp  cos  6  -  j)< 


The  expression  (39)  for  u  .  -  diff  becomes  infinite  on  the  axis  p  ■  0. 
This  is  due  to  the  fact  that  the  axis  is  a  caustic  of  the  singly  diffracted  rays. 
In  order  to  correct  (39)  so  that  it  yields  a  finite  value  on  the  axis,  we  make 
use  of  an  exact  solution  of  the  reduced  wave  equation  which  has  an  axial  caustic. 
This  is  done  in  Appendix  17  ;   it  turns  out  that  (39)  can  be  corrected  on  and  near 
the  axis  by  multiplying  it  by  -s   jfenkp  cos6  sec(kp  cos6  -  y)   J  (kp  cos  6).  J  is 
the  Bessel  function  of  zero  order.  Then  (39)  becomes 


30a  - 


30- 


Figure  ID 

The  far  field  diffraction  pattern  of  a  circular 
hole  of  radius  a  due  to  a  normally  incident  plane 
wave }  ka  =  3n.     The  solid  curve,  based  upon  eq.(38), 
results  from  single  diffraction  with  either  u  -  0  or 
ux  -  0  on  the  screen.     The  dashed  curve,  based  upon 
eq.(52),   results  from  single  and  multiple  diffraction 
for  the  case  u  =■  0  on  the  screen.     The  dotted  portion 
near  0-0  contains  the  caustic  and  shadow  boundary 
corrections  to  both  (38)  and  (52);   it  is  based  upon 
U2;  or  (55) s  which  yield  indistinguishable  curves. 
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ik(x2+a2)^ 


ya  cos  6         [       ,6  ^  ii\  ,6  A  n 

Single  diff  =  *  2(^2}lA   LSSC(?  +  *>  i  CSC(2  +  H 


J  (kp  cos  6)e 


which  is  finite  on  and  near  the  axis. 

Far  from  the  screen,   cos  6    %   a/x  and  then  (bP)  becomes   (with  sin  0  =  p/x) 

2 
(la)       Usingle  diff  =  "  Ie  *MMx  +  |^)]  [sec(|  +  J)  ±  csc(|  +  £ )]jo(ka  sin  0). 

We  cannot  use  (Ul)  directly  to  compute  the  transmission  cross  section  of  the 
aperture  because  6  approaches  n/2  as  x  becomes  large,   and  then  sec  (■*•  +  j-)  ~    2x/a 
becomes  infinite.     This  is  due  to  the  fact  that  the  shadow  boundary  extends  in 

the  direction  6  =  •*,  which  is  parallel  to  the  axis.     Instead  we  first  combine 

iloc 
usingle  diff  with  ugeom»  which  is  e         in  the  lit  region.     If  we  write 

exp|ik(x  +  |^)    «^eikx(l  +  ^-)  then  we  find  from  (hi)  that 

(U2)       u  +u.     .       ..«    «  e1Jac-4=ei3ac(l  +  ^)S£  J  (kasinjf) 

geom        single  diff  2x  2x       a     ox  r' 

■ |jg^-  •  ina    JQ(ka  sin  0). 


From  (1±2)  we  see  that  u     and  u  .  ,   ■,.-.«  combine  near  the  axis  and  far 

geora     single  diff 

from  the  screen  to  give  a  spherical  wave  with  a  finite  coefficient.  Thus 
f  (0)  given  by  (38)  is  to  be  replaced  by  ina  J  (ka  sin  0)  for  0   near  zero. 

According  to  the  cross  section  theorem,  the  transmission  cross  section  is  the 

o 

imaginary  part  of  f  (0),  which  is  noW  just  na  .  This  is  the  geometrical  optics 

s 

cross  section. 

In  order  to  obtain  a  first  correction  to  the  geometrical  optics  cross 
section,  we  wi.11  take  account  of  doubly  diffracted  rays.  Each  incident  ray 
gives  rise  to  a  diffracted  ray  in  the  plane  of  the  screen.  This  ray  is  normal 


-  32  - 


to  the  edge,   and  therefore  it  goes  through  the  center  of  the  aperture  and  is 
normally  incident  on  the  opposite  edge.     The  field  on  this  ray,  when  it  reaches 
the  opposite  edge,  is  given  by,   say,   the  first  terra  in  (36)  with  r,   =  2a,   and 
©1  -  5  .     This  field  is 

2ika  -  inA 

(U3)  -S (/?+  /2)   . 

U  /nEa 

From  (U3)   we  see  that  this  field  is   zero  for  a  screen  on  which  the  normal 
derivative  of  the  field  vanishes  as  was  pointed  out  previously.     Thus  we  will 
consider  for  the  present  only  screens  on  which  the  field  itself  vanishes. 
If  we  treat  (U3)  as  the  incident  field,   (18)  yields 

ik(ri.  +  2a)  61   r  1  -1-1/2 

ihh)         Udouble  diff  =  -777= Se°  "2"   Lrl(l  "  3     rlC°S  VJ 


2nk  /a 


eik(r2  +  2a)  62    r  _n  -,-1/2 

+ sec 


2nk  /a 


2 


r0(l  -  a"  r0cos  60) 


.2vu.  -  =     x2^o  ~2, 


From  (kh)  and  (39)  we  see  that  near  the  axis 

tm                                         2e2ika  -  ±j  j-  ,-1 

(U5)       udouble  diff  =  "  -yjj=j= sec  2    [sec(?  +  tf  +  csc(2  +  H}. 


usingle  diff 


Therefore  u,     ,,      ,.„„  becomes  infinite  on  the  axis  just  like  u   .     n      ,.  „,.  does, 
double  diff  single  diff  * 

because  the  axis  is  also  a  caustic  of  the  doubly  diffracted  rays.     It  can  be 
corrected  near  the  axis  in  the  same  way  as  before.      The  result  is  easily  obtained 
by  using  (hi)   in  (US),  which  yields  for  points  near  the  axis  and  far  from  the 
screen, 
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(U6)       udouble  diff —     I  "  k  /  T  j    e  Jo(ka  sin  *> 


Applying  the  cross  section  theorem  to  (U6),  and  adding  in  the  cross  section 
na     obtained  from  (U2)  we  have 


(U7) 


o-        ,       2sin(2ka  -  ft 


na5         "  /5  (ka)3/2 


_  Q   /O 

The  correction  to  the  geometrical  cross  section  is  of  order  (ka)     '     and  is 
oscillatory.     The  result  (U7)  has  also  been  obtained  by  H.  Levine'-    J. 

A  graph  of  c/na     as  a  function  of  ka  is  shown  in  Figure  11,  based 
upon  (U7).     For  comparison  values  of  o-/na     computed  from  an  exact  solution 
of  the  boundary-value  problem  for  the  reduced  wave  equation  are  also  shown. 

We  will  now  calculate  the  multiply  diffracted  field  u    ,  .    , . „„ .     This 
calculation  is  almost  exactly  the  same  as  that  for  the  slit,   since  each  multiply 
diffracted  ray  is  normal  to  the  edge  which  it  hits.     In  fact  the  previous 
calculation  applies  identically  if  we  replace  the  function  F(s,a,d)  in  (29)  by 

n 


2  /2iScI 


(U8)       F(s,a,Q)  =  -  e  |sec  |(©  -  a)  +  esc  |(9  +  a)|    Kl  -  a_1s  sin  0)1 


The  result  (31)  for  u     ,,    ,.„«  now  applies  with  F  given  by  (U8).     As  we  are 
considering  only  normal  incidence,  a  =  0   and  (31)  simplifies  to 

At  points  far  from  the  aperture,   {h$)  becomes 


_o 


-0> 


00 


-N 


-iO 


o 


-in 


-«■ 


-ro 


-OJ 


-3U- 


(Sb) 


u 


mult  diff 


j  try* 

^—  f  (0). 
2nr   m 


The  diffraction  pattern  f  (0)  due  to  multiple  diffraction  is  given  by 


(5D    fm(0)  =  -u  fss  k"3/2  [2  yjas  e-i2ka  +  in^  *  1] 


-1 


cos 


cos(ka  sin  0  -  -r)  -  i  sin  £  sin(ka  sin  0  -  t-J 


cos  0  ysin  p  . 


The  total  diffraction  pattern  f(0)  =  f  (0)  +  f  (0).     From  (5l)  and  (38)  we  find 

s  m 


r 


(52)        |f (0)|   =       jfcia  k"3/2  (sin  0)"1/2  - 


cos 


COS 


sin 


<u      2v  cos  V  cos  V^  -  2p.  sin  ^  sin  >| 


y^      2u.  cos  £  cos  >J  +2v  sin  £  sin  to 

00s  0 
2)1/2 


sin 


1 


cos 


In  (52)  we  have  used  the  abbreviations  Vj  »  ka  sin  0  -  t-  and 


U  /Sea"  cos(2ka  -  £)+2 

u.  =  '■ 

(53)  Unka+U  /idea  cos(2ka  -  £)+l 

h  /iuca  sin(2ka  -  £) 

v  « 

Unka+U  /nka  cos(2ka  -  ■?)+! 


It  is  to  be  noted  that  f  (0)  becomes  infinite  on  the  axis  0  ■  0,  which 
is  a  caustic  of  the  multiply  diffracted  rays,  and  also  on  the  screen  0  ■  -x   , 
which  is  like  a  shadow  boundary.  To  compute  the  transmission  cross  section  by 


the  cross  section  theorem  we  must  correct  f  (0)  at  0  ■  0.  To  this  end  we  again 

m 
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observe  that  near  the  axis  u  ,.  ..--  becomes  infinite  just  like  u  .  ,   ,  .-j. 

mult  diix  single  daff 

and  u,  .  ,  H'ff*  ^n  ^ac^  by  comparing  (U9)  and  (Ui)  we  see  that 
™   umult  diff  "  t1  -  F(2a,n/2,n/2)]-1  udouble 


Thus  we  can  obtain  u  , ,  diff  near  the  axis  and  far  from  the  screen  by  using 

(U6)  in  (51*).  This  yields  the  following  value  for  f  (0)  near  0=0: 

m 

(55)   fm(0)  -  -  J«{S  [2  /55  e"i2ka  +  in/l  ♦  l]"1  J0(ka  sin  0). 

Now  applying  the  cross  section  theorem  to  {$$)   and  adding  the  contribution  na 
from  (U2),  we  obtain 


(56) 


o- 


2 


2  sin  (2ka  -  f)  r        -1/0  i"!"1 

1 _,0  U   1  +  (nka)"1/2  cos(2ka  -  f)  +  (Unka)"1   . 

v£  rv«v/2    L  u       j 


na        /n  (ka) 


The  cross  section  given  by  (56)  is  shown  as  a  function  of  ka  in  Figure  11, 
along  with  the  result  (1±7)  due  only  to  single  and  double  diffraction. 

We  see  from  (U9)  that  u  ,.  d-ff  is  'infinite'  if  the  factor  in  the 
denominator  becomes  zero.  Just  as  in  the  case  of  the  slit,  this  determines 
eigenvalues  of  k  which  are  complex.  The  equation  they  satisfy  is 
2  /luca  e"l2ka  +  in/^  +  1-0.  None  of  the  roots  of  this  equation  has  a 
positive  imaginary  part,  so  they  all  yield  modes  which  decay  in  time. 


-36  - 


9.  Diffraction  by  a  convex  aperture 

Suppose  that  the  rim  of  the  aperture  is  a  smooth,  closed,  convex  curve. 
Then  for  a  normally  incident  plane  wave  it  is  easy  to  find  the  singly  diffracted 
rays  through  any  point  P.  There  are  just  two  such  rays,  and  they  come  from  the 
two  points  on  the  rim  which  are  nearest  to  and  farthest  from  P.  Thus  the  singly 
diffracted  field  contains  two  terms,  so  it  is  of  the  form  (36).  However  the 
radius  a  in  each  term  of  (36)  now  denotes  the  radius  of  curvature  p  of  the  rim 
at  one  of  the  points  of  diffraction.  In  case  the  rim  has  a  center  of  symmetry, 
as,  for  example  an  ellipse  does,  these  two  radii  are  equal.  The  field  far  from 
the  aperture  can  then  be  simplified  to 


ikr 


<*>   "single  diff <P>  -  "  - 


P(Y) 


2nk  sin  0 


sin[ka(Y)sin  0- j]       cos[ka(Y)sin  0-£] 


sin 


cos 


r 


In  equation  (5U)>  which  applies  to  a  scalar  field  satisfying  the  reduced 
wave  equation,  the  upper  sign  refers  to  a  screen  on  which  u  =  0,  while  the  lower 
sign  refers  to  one  on  which  u  =  0.  The  origin  of  the  polar  coordinates  r,  0  and  T 
of  P  is  the  center  of  the  aperture.  The  angle  0   is  the  angle  between  the  x-axis  and 
the  radius  vector  from  the  origin  to  P.  The  angle  Y  is  the  angle  between  the  y-axis 
and  the  projection  of  this  radius  vector  onto  the  plane  of  the  screen.  The  radius 
of  curvature  of  the  aperture  edge  is  denoted  by  p(l)  and  the  width  of  the  aperture 
in  the  Y  direction  is  denoted  by  2a(Y).  By  the  width  in  the  direction  Y  we  mean 
the  distance  between  the  two  tangents  to  the  edge  which  are  perpendicular  to  the 
Y  direction,  and  p(Y)  is  the  radius  of  curvature  at  either  point  of  tangency. 

Consider  now  an  elliptic  aperture  with  a  major  axis  of  length  2b  parallel 
to  the  y-axis  and  a  minor  axis  of  length  2  e  b.  Such  an  aperture  has  a  center 
of  symmetry  so  ($U)  applies  to  it„   The  curvature  and  width  of  this  aperture 
are  easily  found  to  be 
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a(Y)  -  eb  [^(AlXJ^AaA)-1]"1^ 


(56) 


When  these  expressions  are  used  in  (5U)  an  explicit  formula  for  the  singly 
diffracted  far  field  of  an  elliptic  aperture  is  obtained. 

By  calculations  similar  to  those  made  in  the  cases  of  the  slit  and  the 
circular  aperture,  we  could  determine  the  multiply  diffracted  field  for  any 
aperture.  However  each  order  of  diffraction  would  require  a  separate  and  lengthy 
calculation.  Therefore  we  will  not  attempt  this,  but  instead  we  will  consider 
only  the  doubly  diffracted  field  produced  by  an  incident  plane  wave.  We  will  be 
primarily  interested  in  this  field  far  from  the  aperture  in  the  forward  direction, 
i.e.  in  the  direction  of  the  incident  rays.  The  reason  for  this  interest  is  that 
the  far  field  in  the  forward  direction  determines  the  transmission  cross  section 
of  the  aperture,  according  to  the  cross  section  theorem.  Furthermore,  the 
geometrical  optics  field  and  the  singly  diffracted  field  just  yield  the 
geometrical  cross  section  of  the  aperture,  i.e.  the  projected  area  of  the 
aperture  [8] .  Therefore  the  doubly  diffracted  field  in  the  forward  direction 
will  yield  the  first  correction  to  the  geometrical  cross  section. 

It  is  interesting  to  observe  that  in  general  only  a  few  doubly  diffracted 
rays  go  in  the  forward  direction.  To  show  this,  let  us  denote  by  P  a  point  on 
the  edg  at  which  an  incident  ray  is  diffracted  and  by  Q  the  corresponding  point 
of  secondary  diffraction.  Then  one  of  the  diffracted  rays  produced  at  P  must 
traverse  the  segment  PQ  and  be  diffracted  again  at  Q.  Now  a  doubly  diffracted 
ray  at  Q  which  goes  in  the  forward  direction  is  parallel  to  the  incident  ray 
at  P.  This  is  possible,  according  to  the  law  of  diffraction,  only  if  the  segment 
PQ  makes  the  same  angle  with  the  edge  at  Q  as  it  does  at  P.  Furthermore  the 
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tangent  to  the  edge  at  Q  must  be  parallel  to  the  tangent  at  P.  Thus  only  those 
singly  diffracted  rays  which  satisfy  this  condition  will  produce  doubly  diffracted 
rays  going  in  the  forward  direction. 

To  illustrate  this  result,  let  us  consider  normal  incidence.  In  this 
case  the  singly  diffracted  ray  PQ  must  be  normal  to  the  edge  at  both  P  and  Q. 
For  an  ellipse  there  are  just  two  pairs  of  points  -  the  ends  of  the  principal 
axes  -  having  this  property.  Thus  only  four  doubly  diffracted  rays  will  go  in 
the  forward  direction  when  a  plane  wave  is  normally  incident  on  an  elliptic 
aperture.  It  is  easily  seen  that  there  are  always  at  least  two  such  pairs  of 
points  whenever  the  edge  is  a  smooth  closed  convex  curve.  One  pair  lies  at  the 
ends  of  a  diameter  of  maximum  width  and  the  other  pair  at  the  ends  of  a  diameter 
of  minimum  width.  Thus  for  an  perture  with  any  such  curve  as  rim,  a  normally 
incident  plane  wave  will  produce  at  least  four  doubly  diffracted  rays  going  in 
the  forward  direction.  In  general  there  will  be  only  four  such  rays,  although 
in  some  cases  there  may  be  more.  Thus  for  example  any  aperture  of  constant 
width  will  yield  infinitely  many  -  one  from  each  point  on  the  edge.  We  have 
already  considered  the  circular  aperture,  which  showed  this  behavior.  On  the 
other  hand,  for  non-normal  incidence  the  circle  yields  only  two  doubly  diffracted 
rays  in  the  forward  direction. 

The  field  on  a  ray  diffracted  at  P  and  then  again  at  Q  can  be  computed 
by  applying  (12)  twice.  The  result  is,  for  an  incident  plane  wave  of  unit 
amplitude 


(57) 


u 


double  diff 


Vp  [aft  ♦  p£  d)]-l/2[s(l  ♦  p-Js)]"l/2  e*t*<»  *  d  *  '] 


In  (57)  s  denotes  distance  from  Q,  d  is  the  distance  from  P  to  Q,  p^p  and  p^-, 
are  the  values  of  p1  at  P  and  Q  respectively,  and  Dp  and  D„  are  the  corresponding 
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values  of  D.  The  angle  8  at  P  is  7  since  the  diffracted  ray  lies  in  the  plane 
of  the  screen  and  therefore  at  Q,  a  =  ^.  If  the  doubly  diffracted  ray  is  parallel 
to  the  incident  ray  then  p  at  Q  equals  p  at  P  and  6  at  Q  equals  n  +  a  at  P.  Now 
using  (13)  for  D  with  the  upper  sign,  we  have  when  s  is  large  and  p^Q  is  finite, 


ik[¥(P)  +  d  +  s]  +  ii(l-n)yr - 

t  <R ,  e *    yf piq  pip 

(58)       udoublediff=     — ~2~~~ /< 


nks  sin  p  cos  a 


dCPlp+  d) 


Here  n  is  0,1  or  2  according  as  none,  one  or  both  of  l+dpTp  and  p,Q  are  negative. 

It  is  the  number  of  times  the  ray  goes  through  a  caustic.  The  angles  p  aid  a 
are  evaluated  at  P. 

We  can  easily  compute  pup  by  using  (11)  and  noting  that  cos  p  =  I  •  t 
where  I  is  a  unit  vector  along  the  incident  ray  at  the  edge  and  t  is  the  unit 
vector  tangent  to  the  edge.  Differentiating  this  relation  with  respect  to 

arclength  along  the  edge  yields  -p  sin  p  =  I«t,  since  I  is  constant  for  an 

*    -1  A 

incident  plane  wave.  Now  t  =  p  n  and  I«n  =  -  sin  a  sin  p.  Thus  opsin  p  =  sina  sin  p. 

Furthermore  6,  which  is  the  angle  between  the  diffracted  ray  and  the  normal  n,  is 
just  5  -  p.  Therefore  (ll)  becomes 

ppsin  p 
(59)   P1P  -  -  — 


1  +  sin  a 


In  Appendix  V  it  is  shown  that 

pn  sin  p[pp  sin  p  -  d(l  +  sin  a)] 

(60)  pl     -  -  -2 L. T 

(1  +  sin  o)pp  +  PQ)sin  p  -  d(l  +  sin  a)J 

Ftob  (59)  and  (60),  or  from  (A20),  we  can  see  that  at  least  one  of  the  factors  l+dpT- 

and  p-Q  is  negative,  so  n  ■  1  or  2.  Upon  inserting  (59)  and  (60)  into  (58)»  and 
summing  over  all  doubly  diffracted  rays  in  the  forward  direction,  we  obtain 
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(61)    udouble  diff  =  £ 


Here  B  is  defined  by 


eik[Y(P)+d+s]+lj(:Un)  /pJT 
nks  cos  a  sin  8      *  °|BT 


(62)    B  »  (1  +  sin  a)[(pp  +  pQ)sin  B  -  d(l  +  sin  a)]  . 


From  (61)  we  can  obtain  the  contribution  of  the  doubly  diffracted  rays 

to  the  cross  section  by  using  the  cross  section  theorem.  We  must  first  remove 

ik(r+YQ) 
the  spherical  wave  factor  -ke       /2nr  from  (61)  and  then  take  the  imaginary 

part  of  the  result.  In  the  spherical  wave  factor  Y  denotes  the  phase  of  the 

incident  wave  at  the  origin.  Proceeding  in  this  way  and  adding  in  the  geometrical 

cross  section,  we  obtain 

(63)    <r-  cr     -  ^  y    L_-  sill  (k[y(P)  -  Y  +d+s-rl+e)/^ 

geom   ^2  ^-  cos  a  snip     ^L       o    iJ  if   d|B| 

In  (63)  we  have  introduced  e  -  (l-n£  which  is  -  5  if  B  <  0  and  as  0  if  B  >  0, 
since  n  -  2  if  B  <  0  and  n  -  1  if  B  >  0. 

Let  us  apply  the  result  (63)  to  a  circular  aperture.  In  this  case  only 

points  at  opposite  ends  of  a  diameter  have  parallel  tangents,  and  a  diffracted 

ray  traverses  a  diameter  only  if  p  »t.  It  follows  that  only  two  doubly  diffracted 

rays  reach  any  point,  and  they  come  from  the  ends  of  the  diameter  in  the  plane 

containing  the  point,  the  center  and  an  incident  ray.  Then  if  a  denotes  the 

radius  of  the  circle  we  have  pp  ■  pQ   »  a,  d  =  2a,  Y(P)  -  Y  =  a  sin  a,  s-r  ■  a  sin  a 

and  B  =-2a  sin  a  (1+sin  a).  If  a  denotes  the  positive  angle  between  an  incident 

ray  and  the  normal  to  the  screen,  then  for  one  of  the  doubly  diffracted  rays  we 

must  replace  the  a  in  these  formulas  by  -  a.  Inserting  these  expressions  into  (63) 

yields  for  the  circle 
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(6U)     o-    =  cr 


geom        k2cos  a  ySSTa 


sin[2ka(l-sin  a)]         cos  [2ka(l»sin  a\\ 
Vl-sin  a  /l  ♦  sin  a 


It  i3  to  be  noted  that  (6U)  is  not  valid  for  normal  incidence,  since  sin  a  in 
the  denominator  is  then  zero.  In  that  case  our  previous  result  (U7)  applies. 

For  normal  incidence  the  result  (63)  simplifies  considerably  since  then 
a  =  0,  B  =  ■£,  T(F)  ■  Y  ,  r  =  s  and  B  =  pp+p0-d.   Jith  these  values  inserted, 
(63)  becomes 


(65) 


cr  =»  er 


geom 


7* 


/  dTp? 


p'PQ-d| 


sin(kd+e) 


Each  term  in  (65)   appears  twice,  once  with  incidence  at  P  and  the  other  time 
with  incidence  at  Q.     If  the  aperture  has  a  center  of  symmetry,  then  p^  =  pn. 

If  the  aperture  is  an  ellipse  there  are  four  terms  in  (65).     One  pair  of 


terms  corresponds  to  the  ends  of  each  principal  axis.     Thus  if  the  length  of  the 
major  and  minor  axes  are  2a  and  2b  respectively,  and  if  e  ■     1  -  (— )  denotes  tb* 

eccentricity,  then  (65)  yields 


(66) 


o-  ■  <r 


geom  -  75 


k' 


£? 


2 

sin  2kb  +  5-I±    Cos  2ka 
e 


The  field  on  each  doubly  diffracted  ray  in  (61)  is  proportional  to  k" 

-1/2 
because  each  edge  diffraction  introduces  a  factor  k  '  .  Therefore  the 

—2 
correction  to  the  cross  section  (63)  is  proportional  to  k   because  another 

factor  of  k"  is  introduced  in  defining  the  diffraction  pattern.  However,  if  the 

caustic  of  the  doubly  diffracted  rays  extends  in  the  forward  direction  (i.e.  if 

p..-  is  infinite)  then  (58)  is  not  valid.  In  such  a  case  the  field  in  that 

direction  is  increased  by  a  factor  of  a  power  of  k  which  lies  between  k  and  k   . 


Therefore  in  these  cases  the  cross  section  correction  is  proportional  to  a  power 
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-2     -3/2  -3/2 

of  k  between  k   and  k    .  The  maximum  power,  k    ,  occurs  in  the  result 

for  a  circular  aperture  with  normal  incidence.  The  same  power  occurs  in  the 

result  for  a  slit,  but  for  a  different  reason.  In  the  case  of  a  slit,  which  is 

-1/2  -1 

two  dimensional,  the  factor  k    ,  instead  of  k  ,  is  introduced  in  the  definition 

of  the  diffraction  pattern.  This,  with  the  k~  due  to  diffraction,  accounts  for 

-3/2 
the  k  '     in  the  cross  section  correction  for  a  slit. 

Let  us  now  consider  screens  on  which  the  normal  derivative  of  the  field 

vanishes.  In  Sections  7  and  8  we  found  that  for  such  screens  the  field  was  zero 

on  the  doubly  and  multiply  diffracted  rays.  That  result  is  not  correct.  It  was 

obtained  because  the  assumption  that  the  diffracted  field  is  proportional  to 

the  incident  field  is  not  valid  for  such  screens  when  the  incident  ray  lies  in 

the  plane  of  the  screen.  The  correct  assumption  in  this  case  is  that  the 

diffracted  field  is  proportional  to  the  normal  derivative  of  the  incident  field  at 

the  point  of  diffraction.  Therefore  instead  of  (12)  we  should  have  for  such 

rays  and  screens 


<"■>      <w  d« -»'  #  [-a-  &>]** •tt("*) 


edge  diff 

The  new  diffraction  coefficient  D*  has  been  introduced  in  (12' )•  It  may 
depend  upon  the  nature  of  the  field,  the  directions  of  the  incident  and  diffracted 

rays,  and  upon  k.  From  dimensional  considerations,  D*  must  be  proportional  to 

-3/2 

k  J     .  It  follows  that  for  the  slit  and  the  circle,  the  contribution  to  the  trans- 

-5/2 
mission  cross  section  from  doubly  diffracted  rays  will  be  proportional  to  k    . 

For  a  general  aperture  it  will  be  proportional  to  k.  The  diffraction  coefficient 

D  will  be  deduced  for  a  scalar  field  satisfying  the  reduced  wave  equation  in 

Fl2l 
Part  III  of  this  article1-  -1  and  applied  to  diffraction  by  a  slit. 
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10.  Electromagnetic  diffraction  by  a  convex  aperture 

Let  us  now  consider  a  linearly  polarized  plane  electromagnetic  wave 
normally  incident  upon  a  perfectly  conducting  screen  containing  a  convex 
aperture.  As  we  have  already  seen,  two  rays  pass  through  each  point  P. 
Therefore  the  singly  diffracted  electric  field  ES(P)  is  the  sum  of  two  terms 
of  the  type  given  in  equation  (12).  The  diffraction  coefficient  D  is  given 
by  equation  (A13).  In(A13)  we  must  set  a   =  0  because  the  wave  is  normally 
incident.  Now  if  the  incident  amplitude  is  unity  and  if  the  incident  phase  is 
zero  at  the  screen,  then  ES(P)  is  given  by 


(67) 


ES(P) 


2  /2t3c 


r.^1  -  p~  r-^in  6^1 


-1/2 


e1     ex 

sec  -x-  +   esc  -x- 


sin  9n  sec 


i[3ec  T  " 

A 


CSC  -x- 


9, 


-  CSC 


-sin  1 


-cos  I 


0 


In  (67)  the  terra  written  explicitly  refers  to  one  ray  through  P  and  the 
dot  indicates  a  similar  term  for  the  other  ray.  The  angle  T  is  the  angle  between  the 
negative  incident  electric  field  and  the  normal  to  the  edge  and  p.  denotes  the  radius 

of  curvature  of  the  edge,  both  at  the  point  of  diffraction.  The  components  of 
E3(P)  are  those  in  the  directions  of  the  tangent  t,  the  normal  n  and  the  binormal  b 
at  the  point  of  diffraction.  The  components  in  the  other  terra  are  in  terms  of  a 
basis  referred  to  the  other  point  of  diffraction. 


-  uu  - 


Suppose  that  the  aperture  has  a  center  of  symmetry.  .Also,  let  the  point  P 
be  far  from  the  aperture.  Then  for  such  a  point  P  the  two  points  of  diffraction 
are  practically  symmetric.  Therefore  the  radii  of  curvature  are  equal  at  these 
two  points  while  t„  and  n~  are  the  negatives  of  t.  and  n,  respectively.  In 
addition  &_»  n  +  0,  9„  «£  n  -  0,  r,  ■»  r  -  a(f  )sin  0  and  r2  ^  r  +  a(Y)sin  0. 
Here  0  is  the  angle  between  the  normal  to  the  screen  and  a  ray  from  the  center 
of  symmetry  to  P.  The  distance  between  the  two  tangents  to  the  edge  at  the 
points  of  diffraction  -  the  width  of  the  aperture  -  is  2a(Y).  When  these 
simplifications  are  made  in  (67)  it  becomes 


(68)   ES(P) 


ikr 


nr 


2nk  sin  0 


i  s^  n  ♦ cos  *t 

sin  ^   sin  £• 


0 


-cos 


0 


.sin  n   cos  >l 
~1 


sin 


-sin  0 


sin  £ 


cos 


cos  >i 
cos  fe 


sin  Y 


cos  Y 


7** 


r  /2nk  sin  p 


.  sin  h  .  cos  H 

1tt*  — r 

sxn  <k      cos  & 


sin  Y 


-cos  0 


.  sin  Vi   cos  h 

Sin  T    ens  £• 


COS 


2  J 


cos  Y 


-sin0[i^2^-^i|cosY 
I  sin  I   cos  £ 


Here  Y>  =  ka(Y)sin  0  -  r. 
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In  (68)  the  components  of  ES(P)  are  given  in  terms  of  the  basis  t,  n,  b 
at  either  point  of  diffraction.  This  basis  may  just  as  well  be  located  at  the 
center.  Its  orientation  is  determined  by  the  condition  that  Y  is  the  angle 
between  the  negative  incident  electric  field  and  n.  Furthermore,  T  is  just  the  angle 
between  the  incident  electric  field  and  the  projection  on  the  screen  of  the  ray 
from  the  center  to  P.  This  is  so  because  this  projection  is  practically  parallel 
to  the  normal  at  the  point  of  diffraction.  Therefore  if  the  y  axis  is  parallel  to 
the  incident  electric  field,  then  r,  $   and  Y  are  just  the  polar  coordinates  of  P. 

From  (68)  we  see  that  in  the  H-plane,  Y  =  n/2,  the  electric  field  has  just 
one  non-zero  component  which  is  normal  bo  this  plane.  Quantitatively  this  component 
is  exactly  the  same  as  the  scalar  field  diffracted  by  an  aperture  in  a  screen  on 
which  u  =  0  (see  equation  (5U)).  Qualitatively  it  is  somewhat  similar  to  the 
field  diffracted  by  a  slit  of  width  2a(n/2)  when  the  incident  electric  field  is 
parallel  to  the  edges.  On  the  other  hand  in  the  E-plane,  T  -  0,  the  magnetic 
field  has  just  one  non-zero  component,  which  is  also  normal  to  this  plane. 
Quantitatively  it  is  the  same  as  the  scalar  field  diffracted  by  a  hole  in  a 
screen  on  which  u  «  0  (see  (5U)).  Qualitatively  it  is  also  similar  to  the 
field  diffracted  by  a  slit  of  width  2a(0)  when  the  incident  magnetic  field  is 
parallel  to  the  edges.  These  conclusions  follow  at  once  from  a  consideration 
of  the  location  of  the  diffracted  rays  and  the  orientation  of  the  incident  field 
relative  to  the  edge  at  each  point  of  diffraction.  In  particular,  for  a  circular 
hole,  they  apply  everywhere  and  not  merely  in  the  far  field  because  at  every 
point  P  the  two  points  of  diffraction  are  symmetric  about  the  center. 

It  is  convenient  to  express  E  (P)  in  terms  of  the  unit  vectors  r,  tf  and 

b 

Y  of  the  spherical  coordinate  system.  The  transformation  equations  for  a  point  P 
in  the  n,b  plane  are  t  ■  Y,  n  --r  sin  0  -  0  cos  0  and  b  «■  r  cos  0  -  $   sin  ft. 
3y  using  these  relations  in  (68)  we  obtain 
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ikr 


sin  ^r   cos 


2 


ES(P)  =  0 
r 


These  equations  show  clearly  that  the  diffracted  field  is  transverse  to  the  r 
direction,  which  is  the  direction  of  the  diffracted  rays  far  from  the  aperture. 
Furthermore,  they  show  that  the  components  of  the  incident  field  in  and  perpen- 

A  A       A 

dicular  to  the  T  direction  give  rise,  respectively,  to  the  I  and  y   components 
of  Es. 

11.  Diffraction  by  a  grating  of  slits 

Suppose  the  screen  contains  a  number  of  parallel  slits,  all  of  the  sane 
width,  2a,  and  uniformly  spaced.  Such  an  array  of  slits  is  a  diffraction 
grating.  We  can  apply  our  results  to  determine  the  field  diffracted  by  such 
a  grating.  In  particular  we  can  find  the  positions  and  intensities  of  the 
maxima  in  the  diffraction  pattern.  These  maxima  are  the  spectral  orders  pro- 
duced by  the  grating. 

Let  us  begin  by  assuming  that  the  incident  rays  lie  in  planes  normal  to 
the  edges  of  the  slits.  Then  the  problem  is  two-dimensional,  so  we  will  restrict 
our  attention  to  the  xy- plane.  Suppose  that  there  are  2N+1  slits  and  that  the 
centers  of  adjacent  slits  are  separated  by  the  distance  b.  If  the  center  of  the 
middle  slit  is  placed  at  y  ■  0,  and  if  the  slits  are  numbered  from  t  •  -H  to 
t  ■  N,  then  the  center  of  slit  t  is  at  tb.  It  is  convenient  to  introduce  the 
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polar  coordinates  r.  and  0.  with  the  origin  at  the  center  of  slit  t  and  with 
the  polar  axis  parallel  to  the  x-axis. 

The  singly  diffracted  field  is  the  sua  of  the  fields  on  all  the  rays  singly 
diffracted  from  any  of  the  slit  edges.  To  compute  it  we  first  add  together  the 
fields  on  the  two  rays  from  the  two  edges  of  each  slit.  This  gives  the  singly 
diffracted  field  of  each  slit.  Then  we  add  together  the  singly  diffracted  fields 
of  the  various  slits.  If  we  consider  the  field  far  from  the  grating,  then  it 
suffices  to  add  together  the  far  fields  of  the  various  slits.  For  a  plane  wave 
incident  at  angle  a,  the  far  field  of  slit  t  is  given  by  equation  (22  )  with  r 
and  0  replaced  by  r.  and  0.,  and  with  the  extra  factor  e"    sma  to  account 
for  the  phase  of  the  incident  wave.  Thus  we  have 


(70)       usingle  diff 


t"N    nr~     lk(r.-tb  sin  a)  +i£ 


In  order  to  simplify  (  70  )  we  make  use  of  the  polar  coordinates  r  and  0 
with  their  origin  at  the  center  of  the  middle  slit.  At  points  far  from  the 
grating  compared  to  the  grating  length,  which  is  2Nb  +  2a,  the  angles  0.   are  all 
nearly  equal  to  0.  Furthermore  r  is  approximately  r  -  tb  sin  0.  When  these 
substitutions  are  made  in  (  70  )  we  have  instead 

ikr+ij  t-H 

lsingle  diff  "   ~  /2nr  J  f.(^  £.K 


(71) 

rr-      Ucr+i.5  "in KM  ♦  i-)kb(sin0+sina)] 

-  -  >>r-   •  f  (0) 1 

-Tj  kb(singf  ♦   sinafl 


sin 


The  same  result  applies  when  the  grating  has  an  even  nuaber  of  slits,  say  2N, 
provided  that  the  origin  is  at  the  midpoint  of  the  grating  and  that  N  ♦  j-  is 
replaced  by  N. 
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The  result  (  71  )  shows  that  the  diffraction  pattern  of  the  grating  is  the 

product  of  the  diffraction  pattern  of  a  single  slit  and  the  grating  factor, 

which  is  the  quotient  of  sines*  It  is  clear  that  the  sane  result  with  f  (0) 

s 

replaced  by  f  (?)  would  be  obtained  if  we  took  account  of  rays  multiply  diffracted 

■ 

between  the  two  edges  of  each  slit.  However  a  different  result  would  be  obtained 
if  we  also  included  rays  multiply  diffracted  between  edges  of  different  slits. 


12.  Conclusion 

The  method  which  we  have  applied  in  this  article  to  diffraction  by  an  aper- 
ture is  a  natural  extension  of  geometrical  optics.  It  augments  geometrical  op- 
tics by  introducing  diffracted  rays  and  by  assigning  a  field  to  each  ray.  Both 
of  these  aspects  can  be  worked  out  by  the  use  of  elementary  mathematics  alone. 
Nevertheless  it  leads  directly  to  results  which  have  only  recently  been  obtained  by 
elaborate  mathematical  methods,  and  to  others  which  have  not  yet  been  obtained  in 
any  other  way. 

In  Part  II  the  various  "  physical  optics"  theories  based  upon  the 
Kirchhoff  method  will  be  examined.  Those  theories  require  the  evaluation  of 
integrals.  It  will  be  shown  that  when  those  integrals  are  evaluated  for  small 
wavelengths  they  lead  to  results  which  have  an  immediate  interpretation  in 
terms  of  our  singly  diffracted  rays.  However  the  quantitative  value  of  the 
field  obtained  from  the  integrals  will  be  incorrect,  as  is  well  known.  It 
therefore  seems  that  these  integral  methods  have  oertain  disadvantages  which 
our  method  by-passes.  They  are  harder  to  use,  yield  quantitatively  inferior 
results  and  cannot  take  account  of  multiple  diffraction.  On  the  other  hand, 
they  do  yield  finite  -  but  incorrect  -  field  values  on  caustics  and  shadow 
boundaries. 
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Mathematically  we  can  interpret  the  field  constructed  by  our  method  fron 
the  viewpoint  of  asymptotic  expansions.  Let  us  suppose  that  the  field  u  is 
governed  by  the  reduced  wave  equation  (  ^  +  k  )u-0.  The  field  which  we 


have  constructed  consists  of  a  sum  of  terms,  one  associated  with  each  family 
of  rays.  Each  term  satisfies  the  reduced  wave  equation  up  to  terms  in  k~  , 
as  has  been  shown  elsewhere  *-  -1,  It  is  to  be  expected  that  a  series  of  addi- 
tional  terms  proportional  to  higher  powers  of  k   should  be  added  to  the  ampli- 
tude of  each  of  the  terms  we  have  found.  This  has  been  done  in  [h  ]  for  other 
kinds  of  diffraction  problems.  The  infinite  series  of  terms  associated  with 
each  family  of  rays  should  be  so  chosen  as  to  satisfy  the  reduced  wave  equation 
formally,  but  need  not  necessarily  converge.  The  sum  of  the  various  series 
corresponding  to  the  various  families  of  rayB  should  then  constitute  an  asymp- 
totic expansion  of  the  solution  of  the  reduced  wave  equation. 

The  fact  that  our  expressions  become  infinite  on  certain  surfacee  can  then 

be  interpreted  as  being  due  to  a  change  of  asymptotic  behavior  of  the  solution. 
These  surfaces  may  therefore  be  called  n  Stokes  surfaces"  ,  in  analogy  with 

the  terminology  in  the  theory  of  ordinary  differential  equations.  Shadow 

boundaries  are  examples  of  such  surfaces.  Certain  caustics  are  also  of  this 

type.  Others  correspond  to  a  different  asymptotic  behavior  of  the  solution  on 

the  caustic  itself,  the  fields  on  the  two  sides  of  it  having  the  same  form. 

An  interesting  feature  of  our  theory  is  its  relation  to  what  we  have  called 

canonical  problems.  These  are  special  diffraction  problems  with  relatively 

simple  incident  fields  and  boundary  shapes.  By  means  of  our  method,  certain 

coefficients  obtained  from  these  special  problems  can  be  used  in  the  analysis 

of  many  other  more  complicated  problems.  Thus  a  new  motivation  is  provided  for 

solving  these  special  problems,  since  their  solutions  now  have  a  more  general 

significance. 
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Appendix  I,  Caustics  of  diffracted  rays 

We  wish  to  consider  an  edge  given  by  the  vector  x(s),  from  each  point 
of  which  a  cone  of  diffracted  rays  emerges.  Let  p(s)  denote  the  half -angle 
of  the  cone  with  its  vertex  at  x(s).  The  axis  of  this  cone  is  the  tangent  x(s). 
If  y  denotes  a  variable  point  on  this  cone,  then  y  satisfies  the  equation 

(Al)       7  -  *<3>  .  i(3)  -  cos  p(s) 

|y-x(s)| 

The  envelope  or  caustic  of  the  diffracted  rays  is  also  the  envelope  of  these 
cones,  since  these  cones  contain  the  diffracted  rays.  To  find  this  envelope 
we  differentiate  (Al)  with  respect  to  s  and  obtain 


(^2)     (y  -  x)  .  x  -  i  -  (cos  p)*  |y  -  xl  -  2-UE  .  i  cos  p 

|y  -  x| 


Now  using  (Al)  in(&2)  yields 

(A3)      (y  -  x)»(x  -  x[log  cos  s]*)  «  x  -  cos  0 

If  the  parameter  s  denotes  arclength  along  the  edge  then 
x  =  1,  x  =  p~  n  and  x  =  t.  Here  p,n  and  t  are  the  radius  of  curvature,  the 
unit  normal  and  the  unit  tangent  to  the  edge,  respectively.  Introducing  this 
notation  into  (Al)  and  (A3)  we  obtain  as  the  equations  of  the  caustic 


(AJ4)      (y  -  x)»[n  -  tp(log  cos  p)'l  -  p  sin2p 


(A5)  (y  -  x)-t   -  |y  -  x|  cos  0 
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From  (AU)  and  (A£)  we  can  determine  the  distance  p,  from  the  edge  to 
the  caustic  along  any  particular  ray.  This  distance  p,  is  just  -|y  -  x|  if 
we  take  account  of  the  direction  in  which  p..  is  measured  positively.  Thus  we 
obtain 


p  sin 
(A6)      P-,  -  - 


,_  _  _  . — . 

p{$  sin  p  +  cos  6 


The  angle  6  introduced  in  (A6)  is  defined  by  the  equation 


(A7)      cos  6  -  (y  ~x)     .  n 

|y  -  x| 


It  is  just  the  angle  betwsen  the  diffracted  ray  and  the  normal  to  the  edge. 
As  an  example  of  the  use  of  (AU)  and  (AJ>) ,  let  us  suppose  that  the 
edge  is  a  plane  curve-  and  that  a  plane  wave  is  normally  incident.  Then 
P(s)  ■  -x   so  p  -0.  Equations  (AU)  and  (a£)  then  become  (y  -  x)»n  ■  p 
and  (y  -  x)«t  -  0.  The  first  of  these  equations  shows  that  the  caustic  contains 
the  curve  of  centers  of  curvature  of  the  edge,  i.e.,  the  evolute  of  the  edge. 
The  second  equation  shows  that  the  caustic  consists  of  the  set  of  straight  lines 
through  these  points, each  normal  to  the  plane  of  the  edge.  Thus  the  caustic  is 
a  cylinder  with  generators  normal  to  the  plane  of  the  edge,  the  cross  section 
of  which  is  the  evolute  of  the  edge.  If  the  edge  is  the  edge  of  a  thin  flat 
obstacle,  then  the  caustic  will  appear  as  a  bright  line  in  any  cross  section 
of  its  shadow.  Bright  lines  of  exactly  this  shape  have  been  observed  by  J.  Coulson 
and  G.  G.  Becknell'-  ■*.    The  special  case  of  a  circular  disc  is  well  known. 
The  evolute  in  this  case  is  just  the  center  and  the  caustic  is  the  axis. 
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which  appears  as  a  bright  spot  in  a  cross  section  of  the  shadow. 

In  the  case  of  oblique  incidence  or  of  a  non-planar  edge,  it  may  be 
shown  from  (All)  and  k$)   that  with  increasing  distance  behind  the  obstacle 
the  caustic  approaches  a  cylinder.  The  generators  of  this  cylinder  are 
parallel  to  the  incident  rays  and  its  cross  section  is  the  evolute  of  the 
shadow  cross  section.  Thus  in  any  normal  cross  section  of  the  shadow,  far 
enough  from  the  obstacle,  a  bright  line  should  be  seen.  This  bright  line 
should  coincide  with  the  evolute  of  the  shadow  boundary.  Such  bright  lines 
were  also  observed  by  G.  G.  Becknell  and  J.  Coulson'-  J.    However,  they  did 
not  observe  the  departure  from  this  shape  which  must  occur  nearer  to  the 
obstacle,  as  all  their  observations  were  made  far  behind  it. 


Appendix  II.  The  Diffraction  Coefficient  D 

It  has  been  noted  that  the  diffraction  coefficient  D  depends  upon  the 
nature  of  the  field  and  the  kind  of  screen.  Let  us  therefore  consider  only 
scalar  fields  governed  by  the  reduced  wave  equation,  and  screens  on  which  either 
u  or  its  normal  derivative  vanishes.  For  such  screens  and  fields  the  two  dimen- 
sional problem  of  diffraction  of  a  plane  wave  by  a  screen  in  the  form  of  a  half 
plane  has  been  solved  exactly.  By  a  two  dimensional  problem  we  mean  one  in  which 
each  incident  ray  lies  in  a  plane  normal  to  the  edge.  By  solved  we  mean  that  a 
solution  of  the  reduced  wave  equation  satisfying  the  boundary  and  other  conditions 
has  been  constructed. 

The  solution  for  either  boundary  condition  can  be  extended  to  apply  to 
the  three  dimensional  case  in  which  the  incident  rays  lie  in  planes  making  the 
angle  p  with  the  edge.  To  this  end  it  is  merely  necessary  to  replace  k  in  the 


-  53  - 


jlfg     CO^    ft 

solution  by  k  sin  p  and  then  to  multiply  the  result  by  e  ■  ,  provided 

that  the  z-axis  is  parallel  to  the  edge. 

Once  the  solution  has  been  obtained  it  can  be  expanded  asymptotically 
in  k,   for  k  large.     If  the  amplitude  of  the  incident  wave  is  unity  the  result  is 


ik(is  cos  p  +  r  sin  p)+  in/U   p        ..  t 

(A7)  u  ~  ugeQm  - |sec  2  O  -  a)+  esc  |(©  +  a)j 

2  /2nkr  sin  5 


Here  u  is  the  field  determined  by  the  incident  and  reflected  rays  and  r 
is  radial  distance  from  the  edge.  The  angles  are  shown  in  Figure  6  and  the 
upper  sign  applies  to  the  screen  on  which  u  ■  0.  This  is  to  be  compared  with 


the  solution  given  bv  the  present  method.  In  this  case  there  is  ,just  one 
diffracted  ray  through  eac 
Then  using  (12)  we  obtain 


diffracted  ray  through  each  point.  Also  p  =»  0  so  it  follows  that  pT  -  0, 


ik(z0cos  p+s) 

(A8)      u  =  u     +  2£ 

geom       y^ 

Here  z  denotes  the  point  of  diffraction  on  the  edge.  From  the  law  of 
diffraction  we  see  that  s  =  (z  -  z  )cos  p  +  r  sin  p  ■  r/sin  p.  Upon  inserting 
the  first  of  these  expressions  for  s  into  the  exponent  in  (A8)  and  the  second 
expression  into  the  denominator,  our  result  (A8)  becomes 


(A9)      u  ~  u    + 


De 


ik(z  cos  p  +  r  sin  p) 


geom     /r/sin  P 

The  result  (A9)  is  seen  to  depend  upon  z  and  r  in  exactly  the  same  way  as  doe3 
(A7)  and  to  agree  with  it  precisely  if  D  has  the  value  given  in  equation  (13). 
A  similar  consideration  can  be  applied  to  diffraction  bv  a  wedge  of 
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angle  (2-n)n  instead  of  a  screen,  which  is  a  wedge  of  angle  (2-2)n  =0.  In 
this  way  one  obtains  the  diffraction  coefficient 


(A10)     D  =  - 


e 


WU         [,  -        «^-l-  .    .        n^^-1 


2  /2nk  sin  p 


[/   n      0-a\   I  /    n      9+a+n\ 

(cos cos  )   +  (cos  —  -  cos  ) 

[\n       n/     \n       n/ 


Here  again  the  upper  sign  applies  to  the  screen  on  which  u  =  0  and  the  lower  sign 
to  the  screen  on  which  the  normal  derivative  of  u  is  zero.  The  angles  are  defined 
as  in  Figure  6. 

Appendix  III.  The  electromagnetic  diffraction  matrix 

The  electromagnetic  problem  of  diffraction  of  an  arbitrary  plane  wave 

fid 

by  a  perfectly  conducting  half -plane  has  been  solved  exactly.^  This  means  that 

an  exact  solution  of  ^felxwell,s  equations  satisfying  the  boundary  and  other  con- 
dtions  of  the  problem  has  been  constructed  for  a  plane  wave  travelling  in  any 
direction.  This  solution  has  also  been  expanded  asymptotically  for  large  k. 

This  same  problem  can  be  treated  by  the  present  method.  In  this  case 
there  is  exactly  one  diffracted  ray  through  each  point  not  lying  on  the  edge. 
Furthermore,  p  -  0  and  p~  -  0.  We  then  find  that  the  electric  field  E  is  given  by 

Iks 

(ill)  E(P)  -  E    (P)  +  - D  E.   (Q)  . 

x  x       geonr  '     N/-     incx 

In  this  equation  Q  denotes  the  point  on  the  edge  at  which  the  diffracted  ray  through 

P  is  produced.  The  distance  from  P  to  Q  is  s,  the  geometrical  optics  field  is 

denoted  by  E_„  (P)  and  the  diffraction  matrix  is  denoted  by  D. 
"     geom  " 

We  now  compare  the  expanded  form  of  the  exact  solution  with  our  result. 
We  find  that  the  two  are  identical  if  D  is  given  by 
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(112) 


-e 


inA 


2sinp/2nk" 


sec  ^(9-o)+csc  ^(&-hi) 


e 


9 


cos  -x  +  sin  -x 

-2  cot  p  £ £ 

cos  -x  +  sin  ■» 


cot  p 


8       0       1  1 

cos  *  -  sin  *   sec  »(9-a)-csc  •xiQ+a) 

__   .  _ 
cos  p  +  sin  ■*        cos  a 


cos  6 


cos  a 


sin  6 
cos  a 


sec  |(e-a)-csc  |(0-kc)    0 


sec  |(9-o)-csc  |(9-ki) 


This  form  of  D  applies  if  E  is  expressed  in  terms  of  the  basis  (t,n,b).  Here  t 
denotes  the  unit  vector  tangent  to  the  edge  at  the  point  of  diffraction  Q,  n  is 
the  unit  normal  to  the  edge  (pointing  into  the  aperture)  and  b  is  the  unit  binormal 
to  the  edge,  i.e.,  a  normal  to  the  screen  pointing  toward  the  dark  side.  The 
direction  of  t  is  determined  by  the  condition  that  these  three  vectors  form  a 
right-handed  system,  i.e.,  n  x  b  =  t. 

The  possibility  of  defining  D  in  such  a  way  that  our  result  coincides 
with  the  expansion  of  the  exact  solution  is  a  partial  check  on  our  method.  It 
also  follows  from  the  explicit  form  of  D  that  E  is  transverse  to  the  diffracted 
ray,  as  we  have  assumed.  The  expression  for  D  given  in  (A12)  is  the  desired 
form  of  the  electromagnetic  diffraction  matrix. 

It  is  to  be  noted  that  (A12)  simplifies  when  p  -  n/2,  to  the  following 


(A13) 


-e 


inA 


2  /2nk" 


sec  -j(9-a)+csc  ^(9-ki) 


cos  8 

COS  CL 


sin  8 
cos  a 


|sec  |(9-a)-csc  |(8+a)| 
|  sec  |(8-a)-csc  |(8+c)l 
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This  form  of  D  can  also  be  obtained  by  reducing  the  electromagnetic  problem  to 
two  scalar  problems,  one  for  E,  and  the  other  for  H, ,  the  components  of  E  and  H 
along  the  edge.  Each  of  these  components  satisfies  the  reduced  wave  equation, 
while  on  a  perfect  conductor  E.  and  the  normal  derivative  of  H.  are  zero. 
Therefore  the  scalar  diffraction  coefficients  in  equation  (13)  apply  to  these 
two  components.  If  E  and  E.  are  computed  from  H.  and  if  the  scalar  diffraction 
coefficients  are  used,  the  result  (A13)  isggain  obtained  for  the  electromagnetic 
diffraction  matrix. 


Appendix  IV.  The  field  on  an  axial  caustic 

As  has  been  pointed  out  before,  the  field  given  by  the  present  method 
becomes  infinite  on  a  caustic*  A  finite  value  for  the  field  on  the  caustic,  and 
a  corrected  value  for  the  field  near  the  caustic,  can  be  obtained  by  the  use  of 
a  canonical  solution,  as  was  mentioned  before.  In  this  method  an  exact  solution 
of  the  equation  or  equations  governing  the  field  is  required.  When  this  solution 
is  expanded  asymptotically,  for  large  k,  it  must  correspond  to  a  system  of  rays 
possessing  a  caustic  of  the  desired  type.  A  comparison  of  the  exact  solution  with 
its  expansion  shows  how  the  expansion  must  be  modified  to  remain  finite  on  the 
caustic*  It  is  assumed  that  the  same  modification  applies  to  any  field  with  a 
similar  caustic* 

Let  us  apply  this  method  to  fields  satisfying  the  reduced  wave  equation 
and  having  an  axial  caustic*  We  begin  with  an  exact  solution  of  the  reduced  wave 
equation 

(Alii)      u  -  eikx  sin  6  Jn(kp  cos  6)  cos  n0  • 

In  (Cli)  J  denotes  the  Bessel  ftnotion  of  integral  order  nj  x,  p  and  0  are  cylindrical 
coordinates,  and  6  is  a  constant.  The  asymptotic  expansion  of  (All;)  for  k  large,  ob- 
tained by  inserting  the  expansion  of  J  into  (AlU)  is 
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(A15)    u  ~   C03  n  j    exp[-ik(p  cos  6  -  x  sin  6)  +  i(n  +  |)§l  * 
y2nkp  cos  6  L   L  J 

+  exp  ik(p  cos  6  +  x  sin  5)  -  i(n  +  «■)» 


cos 


n  0    o  ikx  sin  6     IT     _   *   .  X\i\ 

f   ■-   2  e         cos  |kp  cos  6  -  (n  +  •r)-x\  • 

cos  6  *-  -* 


ysi^p" 

The  first  form  of  (Al£)  shows  that  u  is  the  field  on  a  set  of  rays  which 
converge  symmetrically  on  the  x-axis,  pass  through  it,  and  then  diverge.  The 
constant  6  is  the  angle  between  any  ray  and  a  normal  to  the  axis,  so  t  -  6  is 
the  angle  between  each  ray  and  the  axis.  The  field  at  each  point  off  the  axis  is  the 
sum  of  the  fields  on  two  rays,  one  going  toward  the  axis  and  the  other  already 
diverging.  Jbr  such  a  set  of  rays,  the  method  we  have  given  for  constructing  the 
field  yields  exactly  the  expressions  given  in  equation  (Al£).  This  can  be  seen 
from  the  right  side  of  equation  (10)  by  making  p.,  infinite.  This  agreement 
provides  additional  support  for  our  method. 

The  expression  for  u  in  (Al£)  is  infinite  at  p  s  0,  while  actually  u  has 
there  the  finite  value  given  by  (Alii).  To  obtain  the  finite  result  in  (Alii)  from 
the  divergent  result  in  (Al£),  we  must  multiply  the  field  given  in  (Al£)  by 

(Al6)    -x  y2rikp  cos  6  sec  kpcos  6  -  (ri  +  *)i  Jn(kP  cos  6)« 

We  assume  that  any  field  governed  by  the  reduced  wave  equation  and  having  an 
axial  caustic  can  be  corrected  on  and  near  the  axis  by  the  same  factor  (Al6). 
The  value  of  n  to  be  used  in  (Al6)  is  determined  by  the  angular  dependence  of 
the  field  under  consideration. 

It  is  to  be  noted  that  the  phase  of  u  given  by  (Al$)  changes  by  exactly 
-  -x   as  a  ray  passes  through  the  caustic.  The  apparent  change  of  n~  is  exactly 
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cancelled  by  a  corresponding  change  ±n  cos  r0   since  0   changes  by  n.     This  result 


agrees  exactly  with  our  equation  (12)  which  yields  a  phase  change  of  -  ^ 


-1 


whenever  a  ray  passes  through  a  caustic.  In  (12)  one  of  the  factors  s  or  1+pT  s 
under  the  radical  in  the  denominator  changes  sign  on  passing  through  a  caustic, 
yielding  the  phase  change*  The  agreement  of  (Al£)  with  this  result  is  another 
check  on  our  theory. 

Appendix  V.  Caustics  of  doubly  diffracted  rays 

When  a  plane  wave  is  incident  upon  an  aperture,  some  of  the  singly 
diffracted  rays  which  it  produces  cross  the  aperture  and  hit  the  edge  again 
producing  doubly  diffracted  rays.  Consider  a  ray  which  is  initially  diffracted 
at  a  point  P  on  the  edge  and  then  diffracted  again  at  another  point  Q.  We  will 
suppose  that  one  doubly  diffracted  ray  leaves  Q  in  the  forward  direction,  i.e. 
in  the  direction  of  the  incident  ray  at  P.  We  wish  to  compute  p,^,  the  distance 
to  the  caustic  of  the  doubly  diffracted  rays,  along  the  ray  which  leaves  Q  in 
the  forward  direction.  This  distance  is  given  by  equation  (ll)  in  which  all 
quantifies  are  to  be  evaluated  at  Q.  We  have  already  seen  that  p  at  Q  is  equal 
to  p  at  P.,  and  that  the  cosine  of  the  angle  between  the  incident  ray  and  the 
normal  to  the  edge  at  P  is  -sin  a  sin  p.  Since  the  diffracted  ray  at  Q  is 
parallel  to  the  incident  ray  at  P,  and  since  the  normal  at  Q  is  the  negative 
of  the  normal  at  P,  it  follows  that  at  Q,cos  6  «  sin  a  sin  p. 

In  order  to  compute  p,  we  note  that  the  center  of  curvature  on  the  ray 
incident  at  Q  is  at  the  distance  d  +  p,p  from  Q.  Now  cos  p  «  I*t  where  I  is  a 

ur.it  vector  along  the  incident  ray  and  t  is  the  unit  tangent  to  the  edge.  Thus 

•        ... 

-p  sin  p  ■  I»t  +  I't  ■  I«t  +  I»n/pA.  As  we  found  in  studying  spherical  waves, 

% 

•  sin  B 

preceding  equation  (20),  I»t  ■    - — r*y  .     Furthermore,  I«n  ■  -sin  p.     Combining 

P1P a 
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these  results  we  have 


(117)  4  sin  P  -  #4  "  ^Ll 

Pi?3  PQ 


Inserting  (A17)  and  the  value  of  cos  6  obtained  above  into  (11)  yields 


(/18)    -i--  -JL,  _  i*2*EL«f 

PlQ         plP+d       PQsinP 

If  (A18)  is  solved  for  p.~  and  if  p-p  is  inserted  from  (59),  equation  (60)  is 

obtained. 

Equation  (A18)  is  precisely  the  mirror  formula  of  optics,  if  we  define 

the  focal  length  fQ  by 

p0sin  p 
(A19)  fQ  "  -  ifesTa 


Then  (A18)  becomes 


1  1 


(A20)  -JL  -  -±-  +  i 

PlQ       Pi?3       fQ 


By  comparing  (A19)  with  equation  (59)  we  see  that  fQ  is  indeed  a  focal  length. 
For  if  a  plane  wave  hits  the  edge  at  Q,  coming  from  the  direction  towards  which 
the  diffracted  ray  goes,  then  fQ  is  the  distance  to  the  caustic  along  the  diffracted 
ray  from  Q  to  P.  The  distance  P-ip+d  is  clearly  the  'object  distance,1  since  it 
is  the  distance  to  the  center  of  curvature  of  the  incident  rays.  Similarly 
P-|0  is  the  'image*  distance  because  it  is  the  distance  to  the  caustic. 
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